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ABSTRACT 


This report is addressed to some aspects of the analysis of vibration- 
induced attitude errors in strapdown inertial navigators. Three common 
pulse torquing techniques for strapdown gyros are described. The im- 
portant considerations in choosing a rebalance scheme are reviewed and 
relative advantages of the three approaches are discussed. Binary and 
ternary pulse rebalanced gyros are analyzed using describing functions. 
Closed loop sinusoidal transfer characteristics are obtained analytically 
and compared with the results of detailed gyro simulations. The effects 
of compliance between the gyro wheel and gimbal are also illustrated by 
analysis and simulation. 

Compensation of limit cycling binary torqued gyros is treated with the 
goal of providing a better closed loop response. A discussion of the dith- 
ered binary gyro is presented, illustrating limit cycle determination and 
closed loop characteristics. A brief discussion of plans for a compre- 
hensive gyro and system error simulation is provided . The appendices 
provide details in support of the report body and also treat diverse sub- 
jects such as dynamic errors in single- degree-of-freedom pendulous ac- 
celerometers, gyro misalignment caused by angular rates about the out- 
put axis and the recovery of angular rate information from pulse torqued 
gyros. The research discussed in this report includes work presented 
in TASC report TR- 101-3 and constitutes an extension of that effort. 
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INTRODUCTION 


1. 1 VIBRATION- INDUCED ATTITUDE ERRORS 
IN ST RAP DOWN INERTIAL SYSTEMS 

The single- degree-of-freedom gyro, originally perfected for 
use in platform inertial guidance systems, can exhibit large errors when 
adapted to the strapdown application. The major sources of these new 
errors are related to the angular motion environment experienced by the 
gyro when it is rigidly attached to a vehicle. For example, the physical 
inertias of several components combine with angular motion to generate 
torques on the float which are erroneously interpreted as angular motion 
about the sensor input axis. Also, the torquing device used to keep the gyro 
pickoff angle small is prone to certain inaccuracies and can create signi- 
ficant instrument errors, some of which depend on the angular motion 
environment. 

In the presence of linear or angular vibrations about several 
gyro axes many of the motion-induced gyro errors can combine to produce 
constant torques on the float. These are interpreted as a constant angular 
rate about the gyro input axis and ultimately as a drift rate about the cor- 
responding strapdown system axis. Other error torques generated by 
angular and linear vibrations do not create constant errors in gyro output. 
However, as a consequence of the manner in which gyro data is processed 
in a strapdown system many of these effects can cause large attitude drift 
rates in the navigator. Prior work, reported in Refs. 1 and 2, indicates 
that many of the individual motion- induced gyro errors can create strap- 
down system drift rates well in excess of those commonly quoted for 
platform inertial navigators. 
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To permit employment of single- degree- of-freedom gyros in 
strapdown navigators it is essential that the vibration- induced system 
errors be analyzed. A more ambitious goal is the ability to synthesize 
or select the best gyroscope for a given mission- related environment. 
Presently, however, neither analysis nor synthesis can be performed with 
confidence because necessary details concerning the behavior of strapdown 
gyros are difficult to describe mathematically. In particular, the trans- 
mission characteristics relating motion- induced error torques to errors 
in the sensor output are complex because the typical inertial grade strap- 
down gyro is coupled with a sampled, nonlinear rebalance device. Accurate 
analysis of vibration-induced strapdown system errors must await develop- 
ment of mathematical descriptions of strapdown gyro transmission char- 
acteristics that are simple and accurate over the frequency range of 
interest. 


In Refs. 1 and 2 an attempt was made to analyze the performance 
of a limit cycling binary gyro. This report extends that approach to cover 
gyros using ternary and time modulation rebalance techniques. The analy- 
ses of the binary and ternary gyros were verified by simulation and several 
comparisons are presented in Sections 4 and 5. It is concluded that the 
analytic approach provides a good means of predicting gyro loop transmis- 
sion characteristics for sinusoidal inputs over most frequencies of interest. 
Subsequent work will treat gyros employing time modulation torquing and 
the transmission of random inputs through strapdown gyro loops. 


1.2 SUMMARY 

References 1 and 2 treat motion- induced gyro error torques, 
and provide some approximate calculations relating the vibration environ 
ment to strapdown system errors. In these two documents an effort was 
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made to describe analytically the transmission characteristics of a binary, 
limit cycling gyro loop. The work discussed in this report continues that 
analysis. Three basic pulse rebalance techniques are described in Section 
2. Important considerations in selecting a torquing scheme are discussed 
and the relative merits of the three approaches are reviewed. A brief dis- 
cussion treating the propagation of vibration-induced errors through a 
linear gyro loop is presented in Section 3. 

The linear response of a limit cycling binary gyro is shown in 
Section 4, to depend on the presence of the unforced oscillation in the loop. 
When inputs of a size and frequency large enough to extinguish the limit 
cycle occur, the binary gyro loop also exhibits a nonlinear response. How- 
ever, the conditions under which this occurs are unlikely and can be pre- 
dicted analytically. The uncompensated binary gyro loop exhibits a highly 
peaked gain at frequencies below that of the limit cycle. It is demonstrated 
that dynamic compensation in the gyro loop can reduce the resonant re- 
sponse or cause it to occur at much higher frequencies. In Section 5 the 
ternary gyro is shown to have a nonlinear closed loop response, particu- 
larly to inputs of high frequency. Both binary and ternary gyro loops were 
shown to suffer a significant deterioration in transfer characteristics if 
sufficient stiffness is not provided in the rotor bearings. This effect is 
of less consequence to the ternary gyro because its nonlinear nature would 
probably preclude use in measuring inputs in the frequency range af- 
fected. 


Motivation is presented in Section 6 for the construction of a 
comprehensive strapdown gyro simulation which evaluates system attitude 
errors. The appendices provide details related to describing function 
analysis and the gyro simulations used, as well as gimbal suspension design, 
dynamic errors in pendulous accelerometers and recovery of angular rate 
information from pulse torqued gyros. 
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Each major section opens with a summary of the research 
performed, a synopsis of the results obtained and a brief discussion rela- 
ting these results to strapdown gyro design. 
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2 . 


DESCRIPTION OF PULSE TORQUING SCHEMES 


SUMMARY 


Three commonly-used nonlinear rebalance techniques 
for strapdown gyros are treated - binary and ternary 
delta modulation and time modulation. Important char- 
acteristics of pulse rebalance schemes are discussed 
and the relative disadvantages of each technique are 
summarized in tabular form. It is concluded that the 
time modulation approach provides better resolution 
than either delta modulation scheme while exhibiting 
no important disadvantages. 


2. 1 PULSE REBALANCED GYROS 

The single-degree-of-freedom gyro consists of a spinning mass 
supported by a gimbal structure which is free to rotate with respect to the 
instrument case about one axis. The basic elements of this device are 
illustrated in Fig. 2. 1-1. The axis of gimbal rotation is usually called 
the output axis and a device is included to measure the output angle, a Q . 

In most high quality gyros the wheel and gimbal are enclosed in and at- 
tached to a cylindrical structure, or float, and a fluid is inserted between 
the float and case. The fluid serves two purposes: it provides a buoyant 
force to unload the gimbal bearings and it damps angular motion between 
the gyro case and float. In addition to the output axis, two other directions 
are specified with respect to the instrument case. They are the nominal 
direction of the rotor spin axis (rotor angular momentum) and the third 
axis required to define an orthogonal coordinate system. The latter is 
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usually called the gyro input axis. Angular motion about the input axis 
creates a gyroscopic torque on the rotor and gimbal assembly, in the di- 
rection of the output axis. Specifically, a positive angular rate about the 
input axis shown in Fig. 2. 1-1 generates a negative torque about the output 
axis. If damping between the case and float is linear and no other torques 
occur, the output angle is proportional to the integral of angular rate about 
a direction normal to the output and rotor spin axes. 


Spin Reference Axis 



Figure 2. 1-1 Single-Degree-of- Freedom Gyro 


As the float rotates with respect to the gyro case, the sensitive 
axis of the instrument moves with it; if the angular rates measured are to 
be expressed in case- fixed coordinates a transformation must be performed 
between float and case axis systems. Also, if a 0 is permitted to become 
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large, a wide angle transducer is required to measure this angle. In most 
applications the output angle is not allowed to exceed a few milliradians. 
Consequently, a very accurate pickoff is employed to measure a 0 and the 
sensitive axis is assumed to be fixed with respect to the instrument case. 

In inertial navigators employing a stable platform, gyros of this 
type are used to sense angular motion with respect to an inertially fixed 
reference and the platform is rotated in a manner which seeks to null the 
output angle. Consequently, the sensitive axis of the gyro is always in 
close coincidence with the case- fixed input axis and a 0 is an accurate mea- 
sure of angular rotation about that direction. When it is necessary to 
rotate platform inertial navigation systems relative to a fixed coordinate 
frame, a mechanism for applying precise torques about the output axis is 
added to the gyro. The gimbal, rotor and float assembly then acts as a 
torque summer, interpreting the moment applied by the torque generator 
in terms of an angular rate about the input axis. Since the platform is ro- 
tated in a manner which always nulls the output angle, a true angular rate 
results, equal in magnitude and of opposite sign from that represented by 
the torque generator output. 

Gyros used in strapdown navigators are likely to experience 
large angular motions about the case input axis. In order to avoid the com- 
plications connected with a large output angle, a powerful torque generator 
is included in the instrument and the rotor spin axis is kept in close coin- 
cidence with the spin reference axis. By keeping the output angle small, 
variations in a 0 can be attributed to angular motion about the input axis. 
Information concerning angular motion about the input axis is largely con- 
tained in the history of torque required to cage, or restrain, the gyro 
about its output axis. The equivalence between torque applied to the gimbal 
and input angular rate is employed in a platform system to generate motion. 
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The same relation can be used in the strapdown application to imply angular 
rate from torque. Figure 2. 1-2 is a schematic representation of a rebal- 
anced single- degree-of-freedom gyro. 


Rebalance 

Single- Degree-of- Freedom Gyroscope Electronics 



Figure 2. 1-2 Rebalanced Gyro Schematic 


Direct determination of the torque applied is difficult and most 
designers rely on measurements of current flow in the torquer windings 
to indicate the moment produced. Unfortunately, the torque generator does 
not exhibit a gain characteristic which is sufficiently linear, and if a wide 
range of torquer currents is employed extensive calibration is required. 

A common way to avoid problems with torquer nonlinearity to permit ap- 
plication of only two torque levels, equal in magnitude and opposite in sign. 
(A third, or zero level may also be considered. ) in this case, the torque 
generator need only be calibrated and adjusted at two distinct points on its 
gain curve, permitting very accurate control of the moments applied to the 
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gyro output axis. A further refinement to the gyro rebalance mechanism 
is provided by carefully dividing the torquer current into pulses of known 
duration. Each pulse then represents an incremental rotation of the gyro 
about its input axis and the gyro rebalance loop output is well suited for 
calculating system attitude changes in a digital computer. In the ideal 
situation, the rebalance torque occurs as a sequence of pulses of positive 
and negative torque, each pulse having the same weight, and the strapdown 
gyro output consists of a train of positive and negative electronic impulses 
accurately representing incremental rotations about the sensor input axis. 
The single- degree-of-freedom gyro is thus said to be "pulse rebalanced." 

Binary Delta Modulation - The binary delta modulation pulse 
rebalance scheme always provides a rebalance torque. The possibility of 
omitting a torque pulse is purposely excluded. This mechanization can be 
represented by a binary logic element inserted between the torquer and a 
signal representing the gyro output angle (see Fig. 2. 1-3). The torquer 
current, of magnitude D, is supplied with a direction (sign) depending on 
the sign of the output angle. Pulses of torque are created by only allowing 
the logic output to change at fixed intervals, T. A further refinement may 
be employed, permitting torque to be applied only over a known fraction of 
each interval, but this has little bearing on the aspects of pulse rebalance 
techniques under consideration here. Figure 2. 1-4 shows a representative 
train of torque pulses applied to a gyro gimbal when binary pulse torquing 
is employed. The shape of the pulses generated differs from a rectangle 
because of lags in both the torquer and the current source driving it. 

Ternary Delta Modulation - The ternary delta modulation techni- 
que differs from the scheme described above by permitting the absence of 
a torque pulse when the output angle magnitude is small. The mechanization 
can be represented by a relay with an input deadzone inserted between the 
output angle signal and the torquer. Figure 2. 1-5 illustrates the ternary 
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Output (T< 


D 


-D 


Figure 2.1-3 Two Level Logic 



Figure 2. 1-4 Representative Torque Output from 

Delta Modulated Pulse Rebalance Technique 
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detection logic. The deadzone width, 26 , is symmetric about zero and 
defines a range of output angles for which no rebalance torque is applied. 
For low frequency inputs 6 can be related to an angular rotation about the 
input axis of the gyro which will not be detected. However, for input fre- 
quencies above the float bandwidth, the relation is more complex. All 
other aspects of the ternary delta modulation technique are identical to 
those of the binary approach just discussed. 


Output 



Figure 2.1-5 Three Level Logic 


Time Modulation - A third pulse torquing technique, time modu- 
lated torquing, adds the gyro output angle to a periodic waveform and em- 
ploys binary logic on the sum. Figure 2. 1-6 illustrates this approach. 

The binary logic is interrogated at a frequency which is a multiple of the 
sawtooth frequency. In this manner, each period of the added waveform 
is divided into an even integral number of parts. Current driving the torque 
generator can only be switched at the moment the binary logic is interro- 
gated. 
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Sawtooth Wave 


Float Angle 



+ 





Two Level Logic 


Figure 2. 1-6 Time Modulated Torquing 


K the output angle is at null, the torquer provides a pair of 
pulses, one positive and one negative. The torque pulses are of equal 
duration and occur at the frequency of the sawtooth wave. The net effect 
of each pulse pair on the output angle is zero, but the float will oscillate 
at the sawtooth frequency. To avoid detecting this float oscillation the 
signal representing output angle is usually passed through an appropriate 
filter. If the filtered indication of output angle differs from null, the pulse 
pairs appear at the same frequency but the pulses in each pair are of different 
duration, is illustrated in Fig. 2. 1-7. The difference between the length 
of two adjacent pulses is such that the net effect is to drive the output angle 
toward null. The sawtooth amplitude is chosen to ensure that the torquer 
always provides a moment in the same direction at the start of each saw- 
tooth period; the amplitude of the added waveform is larger than any anti- 
cipated signal representation of a 0 . For this reason, the torquer current 
always changes sign at the start of a new sawtooth. Also, the current is 
only allowed to reverse once during each period of the added waveform. 

These restrictions ensure the appearance of one positive and one negative 
pulse, in the same sequence, each sawtooth period. 
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Figure 2„ 1-7 Torquer Waveform for 

Time Modulated Torquing 


When time modulated torquing is employed, the gyro output can 
take two forms. A pulse sequence can be provided at an even multiple of 
the sawtooth frequency, corresponding to the. rate at which the binary logic 
is interrogated. In this case, the output would appear as a long series of 
pulses of one sign followed by a long series of pulses with the opposite 
value. Each pulse represents an incremental rotation about the input axis, 
of a magnitude determined by the torquer capability and the interrogation 
frequency. Most of the information contained in this form of output repre- 
sents the added waveform. Only the difference between the number of pulses 
in two successive series can be construed to represent a 0 and consequently 
motion about the gyro input axis. Commonly, this difference is determined 
by summing pulses over each period of the sawtooth. The net difference, 
available at the frequency of the sawtooth wave, represents the information 
content of the gyro output. Since the sawtooth typically has a one milli- 
second period this alternate form of the gyro output has a bandwidth com- 
patible with many of the computers used to process angular motion infor- 
mation. 
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If a fixed increment Digital Differential Analyzer (DDA) is used to 
process gyro outputs, the high frequency pulse train can be employed. 
However, as stated earlier most of the information contained is related 
to the oscillatory signal that was added to a Q and the attitude calculations 
will exhibit this oscillation. If all gyros are not provided with the same 
sawtooth wave serious pseudo- coning system drift rates can be caused. 

Summary - All three pulse torquing techniques provide digital 
outputs which represent increments of angular rotation about the sensor 
input axis. The angular rotation is implied from a carefully calibrated 
torque pulse applied to the gyro gimbals when a pulse occurs. Any devia- 
tion of the actual time integral of torque from that represented by each 
output pulse will cause an error which is not recoverable, i.e. , differs 
in nature from gyro output errors caused by storage of information by the 
float. Constant deviation of the torque pulse can cause constant gyro drift 
rates while random variations provide random output errors. Differences 
between the net torque pulse generated and that implied by each output 
quantum are a major source of errors in strapdown sensors. 

When system accuracy is considered, factors such as quantiza- 
tion and the information rate of pulse rebalanced gyros are important. 

The presence of unforced oscillations in the gyro output is also considered 
when comparing strapdown gyro torquing schemes. In addition, frequency 
response characteristics and energy dissipation in the gyro loops are sig- 
nificant. The major causes of pulse rebalancing errors are described 
below, followed by a discussion of the relative advantages of each torquing 
scheme. 
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2.2 ASPECTS OF PULSE TORQUING GYRO ERRORS 


Torque Pulse Variations - Frequently the magnitudes of torque 
pulses applied to the float differ from those inferred by the corresponding 
electrical impulses at the gyro output. If the deviations are constant, 
calibration of the strapdown gyro loop will eliminate some of the errors 
generated. From Fig. 2. 2-1 it can be seen that a bias error can result for 
both the binary and ternary pulse rebalanced sensors. Also scale factor 
errors are present in both instruments. However, removing errors of 
this kind from the ternary gyro generally requires two scale factor correc- 
tions. The correction used at a particular time depends on the sign of the 
input angular rate. This represents an unwanted complication in proces- 
sing the gyro output. When the rebalance electronics drift, causing uncom- 
pensated variations in torque pulse weights, new bias and scale factor errors 
in the binary gyro result. The bias error is caused by the fact that torque 
pulses are always being generated in the binary rebalance loop. It is in- 
dependent of the angular motion environment. On the other hand, only 
scale factor errors are generated in the ternary loop when the rebalance 
electronics drift. Consequently, if the input angular rates are small the 
ternary gyro will be less sensitive to unknown variations in torque pulse 
weights. 


Quantization - The output of each pulse rebalanced strapdown 
gyro is a series of digital pulses representing incremental rotations about 
the input axis. This form is not well suited for drift-free calculation of 
the direction cosine matrix because information is lost regarding the order 
in which rotations take place about the system axes. The commutativity 
errors which can result are largely related to the gyro output quantization 
level and the manner and speed with which the outputs are processed. 
(Reference 4 provides an approximate expression for commutation errors 
generated by gyro output quantization. In the relation stated there, the 
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(b) Ternary Torquing 


Figure 2.2-1 Drift Rates Resulting from Incorrect 

Torque Pulse Weights 
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drift rate about a particular axis is given by half a gyro output quantum 
multiplied by the difference in angular rates experienced about the other 
two axes. ) As a general rule, fine quantization is to be preferred. 

Delta modulation pulse torquing schemes achieve fine quantization 
by interrogating the nonlinear logic element at a high frequency, premitting 
the torque level to change sign rapidly. While the size of the output quantum 
in the time modulation technique is also determined by the interrogation 
frequency, the rate at which torquer current changes sign is determined 
by the frequency of the sawtooth wave. In gyros employing delta modula- 
tion rebalance schemes, quantization is limited by time constants in the 
torque generation mechanisms. If the pulses are of too short a duration, 
the torquer current will not reach its design level before the next potential 
switching time occurs. Because of eddy current effects, when the torquer 
lag is on the order of the pulse period, two consecutive torque pulses in 
the same direction will not have the same time integral (effect on a G ) as 
the sum of two separate pulses (see Fig. 2.2-2). Consequently, all cur- 
rent pulses will not have equal weight and an accurate measure of angular 
motion about the gyro input axis cannot be found by counting output pulses. 



Figure 2.2-2 Torque Pulse Variation when the Torquer Lag 

is a Significant Fraction of the Pulse Period 
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In the time modulation torquing scheme, quantization is not tied 
directly to torque pulse width. The current switching frequency is strictly 
limited and quantization is determined by the number of points in each 
cycle of the sawtooth at which the torque can change sign. In the absence 
of switching delays, infinitesimal quantization is possible. The practical 
limit is imposed by the ability to switch the torque levels accurately in 
time. Using this approach, torquer time constant problems are eliminated 
by choosing the amplitude of the added waveform, so that torquer current 
always reaches the desired level before the next change occurs. 

Information Rate - The information rate, defined here as the 
frequency with which new indications of system angular motion will be avail- 
able to the attitude transformation computer, can depend on the pulse 
torquing technique chosen. Typically, a gyro with time modulation rebal- 
ance torquing can provide output impulses faster than one using either of 
the delta modulation approaches. Each impulse is interpreted as an incre- 
ment of angular rotation about the instrument’s input axis. However, as 
discussed earlier, most of the information actually present in the raw 
output of this rebalance technique relates to the oscillatory signal added 
to a 0 and not to angular motion experienced by the gyro. The information 
rate of a gyro employing time modulated torquing is essentially limited by 
the frequency of the additive oscillatory signaLused in mechanizing this 
approach. It would then appear that the information rate of gyros using 
the delta modulation approach to rebalance torquing is higher than those 
employing the time modulation technique, when present interrogation fre- 
quencies are considered. However, the practical rate at which angular 
motion can be sensed by this kind of strapdown gyro is usually limited by 
lags in the float dynamics, not the sampling frequency; output pulses which 
change sign at a rate near that of the interrogation clock are more likely 
to be caused by electrical noise than high frequency motion inputs, hi 
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summary, the useful information rate of a contemporary floated single- 
degree-of-freedom strapdown gyro generally is limited not by the pulse 
rebalance technique employed but by float dynamics. 

Unforced Oscillations - All three pulse torquing schemes pro- 
duce unforced oscillations in the gyro loop at one time or another. The 
limit cycle behavior of a control loop containing a binary nonlinearity is 
well known (Ref. 1). In the presence of many inputs, the three level 
pulse torqued gyro and the time modulation scheme will also exhibit cyclic 
errors. The unwanted oscillations all result from the use of torque pulses 
to rebalance the gyro. Limit cycles in pulse rebalanced gyros can be mis- 
taken for coning motion if they occur with proper phase and at the same 
frequency. System attitude drift rates are generated in this manner. The 
binary delta modulation scheme is usually thought to be the least satisfac- 
tory from this point of view. However, it can be shown that slight differ- 
ences between limit cycle periods of the different gyros in a strapdown triad 
will prevent large system errors of this type (see Ref. 1). Thus, gyros 
with slightly different characteristics may be desired. 

Gyro Frequency Response Characteristics - Closed loop response 
to sinusoidal inputs is a useful description of behavior. Closed loop gain 
and phase shift characteristics are useful in determining strapdown system 
errors when the gyro loop is assumed to be linear (see Section 3). Approxi- 
mate determination of these characteristics is discussed in this report for 
all three rebalance techniques. Torque loops employing binary logic ap- 
pear to be more readily described in terms of linear behavior than those 
using ternary logic (see Section 5. 1), offering an advantage when analyzing 
and compensating dynamic system errors. 
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Energy Dissipation - Temperature gradients within the gyro 
structure are a major cause of single-degree-of-freedom gyro drift. K 
the gradients can be held constant, proper testing will enable accurate 
compensation of these effects. While it consumes less power in most mo- 
tion environments, the ternary delta modulation torquing scheme does not 
usually provide a uniform energy flow into the torquer. To avoid this 
obvious cause of varying temperature gradients, the current source can 
be switched into a dummy load when no torque to called for. However, if 
the same heat pattern is to be produced, this requires a dummy winding 
the same size as the torquer. Situations do exist where the ternary tech- 
nique will provide essentially constant energy flow without the extra wind- 
ing. For example, while operating in the angular rate environment 
generated by a limit cycling spacecraft attitude control system. 


2.3 COMPARISON OF PULSE TORQUING TECHNIQUES 

Binary Delta Modulated Torquing - The binary delta modulation 
rebalance technique provides a gyro loop response which is linear because 
of the limit cycle always present in this kind of system. However, some 
compensation within the gyro loop may be desired to improve its behavior. 
Because torque is always being generated, energy is dissipated within 
this kind of an instrument at a constant rate and no significant errors are 
caused by variations in temperature gradients. A high rate of power 
consumption results, even in the most benign motion environment. 

A limit cycle is always present in this instrument and consequently strap- 
down system attitude drift rates can result. The output quantization level 
of a gyro using this rebalance technique is limited by torquer dynamics 
and the output resolution is generally poorer than can be obtained from 
sensors using time modulation torquing. Also, the instrument is relatively 
sensitive to uncompensated deviations in the size of rebalance torque pulses. 
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Binary delta modulation requires a minimum of electronic components, 
suggesting better reliability. Section 4. 1, 4.2 and 4. 3 of this report treat 
aspects of binary pulse torquing. 

Ternary Delta Modulated Torquing - Gyros employing ternary 
delta modulated torquing can conserve electrical energy in a benign envir- 
onment such as that encountered during long space missions. In addition, 
there is less tendency to exhibit cyclic output errors. Because a null 
torquer current can be achieved with high accuracy, they usually have a 
lower drift rate in the absence of an input and are less sensitive to devia- 
tions in torque pulse size. However, in a variable environment unstation- 
ary temperature gradients can be generated inside the instrument unless a 
dummy electrical load is used to dissipate energy when no torque pulse is 
called for. If this additional device is added to the gyro, low power con- 
sumption is sacrificed and the weight and size of the instrument are in- 
creased. As with the gyro using binary delta modulated torquing, the 
ternary rebalanced sensor suffers from an output resolution that is limited 
by torquer dynamics. A gyro loop using this torquing scheme exhibits a 
nonlinear response to high frequency inputs, casting doubt on its value as 
a strapdown sensor in vigorous angular vibration environments. The 
ternary gyro also displays multi-valued response characteristics in the 
high frequency range. Section 5 of this report deals with the ternary pulse 
rebalanced gyro. 

Time Modulated Torquing - The output resolution of a gyro 
employing time modulated torquing is essentially independent of torquer 
dynamics. Consequently, this type of instrument uses finer quantization 
than is commonly found in other pulse rebalanced gyros. An oscillatory 
error can appear in the gyro output but this is much smaller than those 
occurring with binary delta modulation techniques. A high power consump- 
tion is characteristic of this rebalance scheme but temperature gradients 
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within the device are held constant. This type of instrument exhibits a 
relatively large sensitivity to variations in torque pulse weight. Analysis 
indicates that a more satisfactory linear response can be obtained from 
this type of gyro than is possible using simple binary torquing in an uncom- 
pensated loop. Of the three techniques discussed, time modulated torquing 
requires the largest number of electronic components. 

Summary - The three pulse torquing techniques discussed are 
compared in Table 2.3-1. ’Some representative values of important para- 
meters are listed to aid in the comparison. None of the approaches dis- 
cussed offers a clear advantage in all situations. If constant energy dissi- 
pation over a wide range of inputs and linear response are required at the 
outset, the ternary approach can probably be ruled out. However, when 
the angular motion environment is quiet and reasonably constant the ter- 
nary gyro has the advantage of low power consumption and is less affected 
by torque pulse deviations. All three techniques appear capable of providing 
the same information rate, but the time modulation approach can give better 
output resolution because it permits output quantization to be specified 
independent of the torquer design. 


22 



TABLE 2.3-1 


COMPARISON OF 
PULSE TORQUE TECHNIQUES 


UNDESIRABLE 

FEATURES 

REBALANCE TECHNIQUES 

COMMENTS 

BINARY 

DELTA 

MODULATION 

TERNARY 

DELTA 

MODULATION 

TIME 

MODULATION 

Output Quantum 
Level Limited by 
Torquer Dynamics. 

X 

X 


Gyro output quantum, A0, is 
independent of torque lags when 
the time modulation technique 
is employed. 

A9 min > 3a Vax T tg. 

about 20 sec or more when 
“imax = 1 rad/sec 

A8< I sec for the 
Hamilton Standard 
RI-1170 under simi- 
lar circumstances. 

Unforced Oscillations 
in the Sensor Output. 

X 


X 

Ternary gyro can also exhibit 
output oscillations. The finer 
quantization available with the 
time modulation schemes 
reduces the amplitude of these 
errors. 

Amplitude Dependent 

Response 

Characteristics. 


X 


See Fig. 5.1-7 for an illustra- 
tion of the amplitude dependent 
response characteristics of the 
ternary gyro. 

Nonstationary 
Temperature 
Gradients in the Gyro. 


X 


This can be avoided in the 
ternary gyro also but at the ex- 
pense of significant increases in 
torquer size and weight. 

High Power 
Consumption in a 
Benign Environment. 

X 


X 

The ternary gyro also suffers 
this problem if a dummy elec- 
trical load is added to give 
stationary temperature gradients. 
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3. 


STRAPDOWN SYSTEM ATTITUDE ERRORS 


ARISING FROM VIBRATIONAL MOTION 


SUMMARY 


The system attitude drift rate caused by the interaction 
of vehicle vibrations and single- degree-of-freedom 
(SDF) gyro parameters is discussed. Major sources 
of system drift rate are developed followed by a dis- 
cussion of possible gyro parameter selections to re- 
duce attitude error. It is assumed that the gyros can 
be represented by a linear transfer function relating 
torque on the float to sensor output. No errors are 
presumed to be generated in the calculations used to 
process gyro outputs. The importance of simple ana- 
lytic descriptions of gyro loop transfer characteristics 
is discussed. 


3. 1 DESCRIPTION OF GYRO AND VIBRATION CHARACTERISTICS 

The error torques on a single- degree-of-freedom gyro are de- 
veloped in Ref. 1. In the discussion which follows these torques are 
assumed to propagate through a linear gyro loop whose Laplace transform 
transfer function is represented as G(s). At any particular frequency of 
input the in-phase and quadrature components of the gyro output are speci- 
fied by decomposing G(s) into real and imaginary parts: 

G(s)| = ReG(v) + j ImG(v) (3.1-1) 

s = ]v 


25 



The correlation characteristics of angular and linear vibrations 
are described in terms of spectral density functions when only one axis is 
considered and the components of cross spectral density, cospectrum and 
quadspectrum , when two axes are considered. For example, when angular 
vibrations are described the spectral density is represented by the symbol 
4>rr( u )> which is the mean product of angular rate oscillations at frequency 
v, about the r axis. When angular motion about two different axes is under 
consideration the cospectral density function Q rs (v) represents the mean 
product of two angular rate oscillations at frequency v, one about the r axis, 
the other about the s axis. Only those components of the oscillations which 
are in phase provide contributions to 4> rs . To describe vibrations that are 
correlated but in quadrature, a quadspectral density function 4>J. S (v) is 
defined. It is the mean product of the angular rate oscillation about the r 
axis and the vibration about the s axis phase shifted tt/2 radians. To illus- 
trate, suppose the spin and input axes of a particular gyro experience the 
following angular rates: 

Wg = a sin vt 

= b sin vt + c cos vt 

Then; 

a a 2 

$ (v) £ fL. 

v ss w 2 

& /,.\ a ab 

WisW) = -g- 

JL< /, \ A 3-0 

hs w) = x 
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Similar spectral density functions which describe the correlation charac- 
teristics of linear vibrations are designated ^ rr (i^), ^ rs (^) and ^’ rs (y). 
It is assumed that all necessary density functions are known. 
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Figure 3.1-1 Orientation of Gyro and System Axes 


In order to relate system error growth to vibration-induced gyro 
errors the orientation of each gyro with respect to the vibration axes must 
be specified. The most general orientation is used here. The strapdown 
system (vehicle body) axes in which the vibration spectral density informa- 
tion is stated are designated x, y and z. The gyro which measures angular 
motion about the x axis has, of course, its input axis in the x direction. 
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The spin and output axes of this sensor are parallel to the y and z system 
axes. The gyro measuring angular motion about the y axis has its spin 
and output axes parallel to the z and x directions, etc. Figure 3. 1-1 illus- 
trates the relative orientation of gyro and system axes. 


3.2 DRIFT RATE CONTRIBUTIONS 

Only constant system attitude drift rates are considered here. 
The causes discussed below represent the major contributors to vibration- 
induced system drift rate. However, the list is not necessarily complete. 


3.2.1 Spin-Input Crosscoupling 

This error is caused by coupling between angular motion about 
the nominal rotor spin axis and a float angle generated by angular motion 
about the gyro input axis. It can result in a constant drift rate component 
in the individual gyro. The float angle is given by 


«o(s) 


-Hwi(s) 
y if G(s) 
Ksg*Hg 


where 

H = gyro angular momentum 
ooj = angular rate about the input axis 
K S g = signal generator gain 
Kjg = torquer gain 


28 



The average value of the vibration-induced gyro drift rate o> s a 0 is given by 


H 


= k~ - k T~ Ws(s) G(s) 

sg tg 


(3.2-1) 


H 


K sg K tg o> 


/ dv [Re G(v) Q is (y) - Im G (v) *’ is (u)] 


3.2.2 Spin-Output Crosscoupling 

This error is caused by angular motion about the nominal rotor 
spin axis and a float angle generated by angular acceleration about the gyro 
output axis. It can result in a constant drift rate in the individual gyro. 

The float angle is given by 


«o 


- r oo “o(s) 
KsgKtg 


G(s) 


where 


loo = moment of inertia of the gimbal 
and rotor about the output axis 


The average value of the vibration-induced drift rate is given by 


-a 0 a> s 


‘oo 


KsgKtg 


W 0 (s) io s ( s) G(s) 


-I 


oo 


K 


S| 


g K tg 


(3.2-2) 


y * 00 vdv [ReG(^) % B (v) + Im G(y) % s ( v )] 
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3.2.3 Anisoinertia 


This error is caused by a mismatch of moments of inertia about 
the spin and input axes, in combination with angular motion about these two 
axes. Spin axis angular motion of the gyro case and gimbal is transferred 
to the rotor through the wheel speed control loop. Consequently, the spin 
axis inertia of the rotating member can be considered coupled to the float 
only at vibration frequencies below the bandwidth of this control loop; the 
anisoinertia is a function of frequency. Anisoinertia error torques can 
cause constant vibration-induced drift rates in the individual gyro. The 
drift rate is given by 

anisoinertia drift rate = - (Iy - Issg) ^i^s 

*ss r 

+ ~ g ~~ Q(s) 


where 


l ii 


l ss 


g 


rotor and gimbal moment of inertia 
about the input axis 

gimbal moment of inertia about the 
spin axis 


I S s = rotor moment of inertia about the 
spin axis 


Rotor inertia coupling about the spin axis is represented by the transfer 
function 


Q(s) 


CO, 


s a + cc r ‘ 
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where 

o) r = wheel speed control loop bandwidth 
The expression for the average anisoinertia drift rate is 

av e ra^anisome rtia j * (I .. . Iggg . ^ Q(ll)) 9 . g(v) (3 . 2 . 3 ) 


3.2.4 Torquer Scale Factor Asymmetry 

This error is caused by a lack of symmetry in the gain (scale 
factor) relating gyro output to the rebalance torque applied. In a gyro with 
linear rebalance electronics , a sinusoidal angular rate about the input axis 

Wj = W sin ut 

will cause a mean drift rate of - eW/ ir if the frequency v is within the loop 
bandwidth and the torquer gain for applying positive rebalance moments is 
(1 + c) times that for negative moments. The same drift rate occurs in a 
pulse rebalanced gyro when positive and negative pulses are similarly un- 
equal (see Ref. 1). For the linear gyro the mean drift rate caused by 
scale factor asymmetry and angular vibrations is given by 

mean drift rate caused by cf /*°° , ,2 l*/ 2 

scale factor asymmetry ~ ^ I J ^ *ii^) (3.2-4) 
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3.2.5 Pseudo- Coning 


The cause of this system level error is the combination of an 
angular vibration correctly measured by one gyro with a vibration- induced 
error from a second gyro. The strapdown navigator erroneously believes 
that it is experiencing coning motion and a system drift rate occurs about 
the third system axis. The major cause of pseudo- coning system errors 
is the sensitivity of the single-degree-of-freedom gyro to angular accel- 
eration about its output axis. This effect provides an error signal of the 
proper phase to combine with the true motion measured by another gyro, 
giving the erroneous indication of system coning motion. The average 
pseudo- coning error per gyro pair (of which there are three in each system) 
is 


PSeU driftr^ - tA* ooO'HGMl 2 (3.2-5) 

o 


3.2.6 Undetected System Coning Motion 

If the angular vibration environment contains frequencies beyond 
the gyro bandwidth, the system can experience coning motion which is not 
detected by the gyro. The mean drift rate due to undetected angular rota- 
tion about the system z axis is given by 


undetected coning 
drift rate 


oo 

-J |^(1- |G(K)|) 2 ^ ty (K) 


(3.2-6) 


Similar expressions can be written for drift about the x and y axes by 
proper permutation of the subscripts on the quadspectral density function. 
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3. 2. 7 Linear Vibration Errors 


Linear vibrations can cause error torques about the output axis 
of an SDF gyro. Compliance in the gimbal and rotor structure is the mech- 
anism for rectifying the effects of linear vibrations to generate constant 
drift rates in the individual gyros. If the structural compliance coefficients 
Kj[ j , relating rotor center of mass displacement in the i direction to accel- 
eration in the j direction, are constant the mean drift rate per gyro is 


linear vibration 
drift rate 



tfoiW + K si ^iiW 


(3.2-7) 


+ 0K SS “ ^ii) ^is^) “ ^io ^osO') 
" K is ^ss(^)l 


where 


m = rotor mass 


Unfortunately, the compliance coefficients can be frequency dependent and 
Eq. (3.2-7) is a very simple representation of linear vibration errors. 
For a more detailed treatment see Ref. 7. 


3.3 SYSTEM DRIFT RATE 

The system drift rate about any axis can be found by referring 
to the gyro orientation and properly summing appropriate drift rate terms 
from Section 3. 2. For example, 
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mean x-axis system drift rate = 


- ^ f dt- |GM| 2 * 


f di/jlte G (u) 


xx 


(*) 


1 1/2 


+ Im G(^) 


«W (*) - lr 4L- ^(V) 
K sg K tg y K sg K tg y . 


w 


H 


[Kgg^tg 


kr*y^ ) + K 


sg^tg 


+ jj (^ii " Issg " ^ss r Q( l/ ))®xy( l/ ) 


(3.3-1) 


*zzM 


|GM| 2 -i(l- |GM|) 2 *^ z M 


m“ 

H 


K so ’J'xzM + K si *xxM + (K U -K SS ) t xy (y) 
+ K lo '■'yzM - K is 'J'yytt') 


Some observations regarding gyro design can be made from Eq. 
(3. 3-1). If # X y ( v ) is significant only in a narrow range of frequencies, 
adjusting the wheel speed regulator to provide a constant value of Q(y) 
over that band will permit elimination of the anisoinertia drift rate 
by proper selection of gimbal moments of inertia. All other errors 
induced by angular vibration are in some way affected by the gyro loop 
response, represented by the expression G(v). Torquer asymmetry and 
pseudo-coning errors can be diminished by reducing the gyro bandwidth. 
However, at least for the linear gyro, reducing bandwidth will increase 
crosscoupling errors and will cause information to be lost (for example, 
there will be more undetected coning motion). Evidently, determination 
of the optimum gyro parameters is not a trivial problem and is intimately 
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involved with the model of gyro transfer characteristics. While specifying 
G(v) is easy when a linear rebalance technique is used, most strapdown 
gyros employ pulse torquing. The nonlinear elements in a pulse rebalanced 
gyro complicate the description of transfer characteristics for this type of 
instrument. Consequently, Eq. (3.3-1) cannot be used in a straightforward 
manner to compute vibration-induced errors or as the basis for selecting 
gyro parameters. The following two sections of this report outline attempts 
at describing the response characteristics of pulse rebalanced gyros in a 
concise manner that will facilitate calculation of vibration-induced errors. 
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4 . 


BINARY PULSE REBALANCED GYRO 


SUMMARY 


Research Performed - An analytic technique is used to 
predict the behavior of binary pulse torqued gyros. In 
particular, closed loop frequency response character- 
istics are determined for limit cycling delta modulated 
gyros and for instruments employing time modulated 
torquing. Most of the effort is devoted to analysis of 
sensors employing the binary delta modulation rebalance 
technique. For these gyros the effects of interrogation 
period, float time constant and rotor bearing compliance 
are investigated and a gyro simulation is employed to 
verify the behavior predicted analytically. The influ- 
ence of various compensation networks on closed loop 
behavior is also discussed. 

In this section the effect of sampling in the rebalance 
loop is represented as a linear, first order lag in the 
gyro loop. Rotor bearing compliance is modeled in 
terms of a linear, lightly damped spring-mass system. 

Results - Because of the limit cycle usually exhibited, 
the binary gyro is shown to have an essentially linear 
response to sinusoidal inputs. Good agreement is 
observed between the closed loop gain and phase shift 
characteristics determined analytically and those 
observed during simulation. Analytic techniques are 
also shown to accurately forecast limit cycle quenching - 
the situation where the input signal extinguishes the 
limit cycle in a binary gyro. It is demonstrated that 
the effect of sampling on gyro loop behavior can be 
studied in terms of an equivalent first order dynamic 
lag. An analysis is also performed for a strapdown 
gyro loop employing a time modulation torquing 
technique. It is concluded that this rebalance 
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scheme can provide more flexibility in achieving a de- 
sirable closed-loop gain characteristic. The analysis 
is not checked by a simulation. 

Conclusions Related to Gyro Design - A tool for ana- 
lyzing binary strapdown gyro designs as demonstrated. 
However, at present it has only been validated for the 
case of sinusoidal inputs. Further work is required 
before the same analytic approach can be used with con- 
fidence to compute vibration-induced strapdown system 
errors. This technique is not limited to a particular 
set of linear dynamics within the loop and accurately 
predicts what range of input frequencies will be dis- 
torted by the gyro. An uncompensated limit cycling 
binary loop exhibits a frequency response that can be 
closely approximated by a linear, second order transfer 
function whose damping ratio and natural frequency are 
given by Eqs. (4. 1-1) and (4. 1-2). Unfortunately, the 
damping ratio of the approximate system never exceeds 
one-half, giving high peaks in the closed-loop gain. It 
is demonstrated that the response characteristics of a 
limit cycling binary gyro can be changed by inserting 
linear dynamic compensation between the signal genera- 
tor and the rebalance logic. 

Rotor bearing compliance can have a sizeable effect on 
the closed loop transfer characteristics of a single- 
degree-of-freedom gyro. Stiff rotor bearings must be 
provided to enable a strapdown gyro to realize its full 
potential. Sampling in the gyro rebalance loop also 
lowers the bandwidth and the highest sampling rate 
compatible with other constraints should be used. 

The float time constant affects binary strapdown gyro 
bandwidth. Generally, reducing this parameter raises 
the natural frequency and damping ratio of the approxi- 
mate loop response model. 

Quenching of the binary gyro limit cycle and, consequently, 
destruction of the linear response characteristics of this 
sensor can be accurately predicted. Since prediction is 
possible, a designer should be able to avoid quenching 
through proper choice of gyro parameters. 
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Finally, the time modulated rebalance technique offers 
greater flexibility in adjusting gyro loop transfer char- 
acteristics. 
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4. 1 ANALYSIS OF THE LIMIT CYCLING BINARY GYRO 

This section discusses analytic techniques for deriving the 
closed loop transfer characteristics of the binary pulse torqued gyro loop. 
Gyros with second order dynamics and those with significant compliance 
between rotor and gimbal are considered in separate sections. Closed 
loop gain characteristics predicted analytically are displayed to illustrate 
the effects of different instrument parameters. 


Rotor & Gimbal 


\ 

1 /c 


K Sg 


D 

— 


) ► 

S(TjS -*■ 1) 



— 






iKtgTj 

it 


A# 



K tg 


V 71 


1 - e" sT 


s 



Torque r Zero Order 

Dynamics Hold 


H = Rotor Angular Momentum 
C = Float Damping Coefficient 
itg = Torquer Current 

Figure 4. 1-1 Binary-Torqued Gyro Loop with 

Second Order Gyro Dynamics 


4.1.1 Simple Gyro Dynamics 

Figure 4.1-1 illustrates a simple representation of a binary 
pulse rebalanced gyro loop. The rotor and gimbal are assumed to behave 
as a single rigid body for angular motion about the gyro output axis. When 
sampling is ignored, the response of the gyro loop illustrated to oscillatory 
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inputs closely approximates that of a linear second order system with 
damping ratio and natural frequency given by (Section 2.4, Ref. 1): 


C 


to 


n 



(4.1-1) 


(4.1-2) 


The linearization of the binary gyro loop is a consequence of the limit cycle 
usually found in a system of this type. Under certain circumstances (dis- 
cussed in Appendix D) an oscillatory input to the gyro loop can cause the 
limit cycle to vanish. In that case, the closed loop characteristics of the 
binary torqued instrument differ significantly from those presented above. 
When the limit cycle is not quenched, Eqs. (4. 1-1) and (4. 1-2) can be used 
to study the effect of varying gyro parameters. Figure 4. 1-2 illustrates 
the behavior of the closed loop gain as the float time constant is varied. A 
similar comparison could be made of phase shift. Changes in Tj are easily 
made by varying the float damping coefficient, C. 

Effects of Sampling in the Gyro Loop - In the above discussion 
the effect of sampling in the loop has been ignored for the sake of simpli- 
city. The sample and hold operation contributes an additional frequency- 
sensitive phase lag to the open loop dynamics. This can cause the limit 
cycle frequency and amplitude to be different from those experienced in a 
continuous binary gyro loop. Consequently, there is a different value for 
the nonlinearity gain to the signal and the closed loop frequency response 
is not the same. 
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Frequency (Hz) 

Figure 4 0 1-2 Closed Loop Frequency Response as a 

Function of Float Time Constant 


The sample and hold transfer function is 


S(j^) 


. vT 
-jvT/2 Sm X 
e vT 

T 


(4.1-3) 


where T is the sampler period in seconds. K UT/2 « 1 because the 
sampler frequency is much higher than other frequencies encountered in 
the loop, the transfer function is approximately 


S(jv) ~ 


1 


1 + j 


vT 

~T 


(4. 1-4) 
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The equivalent closed loop natural frequency and damping ratio given in 
Eqs. (4. 1-1) and (4. 1-2) can be derived including the effect of the sampler 
lag expressed by Eq. (4. 1-4). The resulting relations are given below. 


co = 
n 



T. + T, , 

i tg 

2r„ 


,T t T tg+f W 

( T i T t g + ? <vw) 2 


(4. 1-5) 
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(4. 1-6) 


The effect of a finite sampling rate on the closed loop gain is illustrated 
in Fig. 4. 1-3 where response curves are shown for a very high sampling 
frequency and a 10 kc rate. It can be seen that sampling in the gyro loop 
tends to reduce the useful bandwidth of the instrument. 
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4.1.2 Gyro Dynamics with Rotor Bearing Compliance 



Figure 4. 1-4 Binary- To rqued Gyro Loop when Rotor 

Bearing Compliance is Considered 


The binary-to rqued gyro loop with significant rotor-to-gimbal 
compliance is in many respects similar to that discussed in Section 4. 1. 1. 
The difference lies in the transfer functions relating torque on the float to 
input angular rate, and float angle to torque (see Fig. 4. 1-4). When linear 
rotor bearing compliance is present, the transfer functions become 


Mj h / 4p\ 2 s 2 + 4pt^s + >4 

“* \s 4 + 4pii 1 s 3 + (2i^ (l + 2 p2)+ w g 2 )s 2 + 4pt$3 + i4, 


(4.1-7) 
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where 
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These relations are derived in Appendix I„ If the coupling between the 
wheel and gimbal is represented by a linear second order system, p and 
v n are the damping ratio and natural frequency of this system. They 
are computed from the stiffness between the rotor and gimbal, k, a 
damping coefficient, d, and the rotor moment of inertia about input and 
output axes, I r . 




( 4 . 1 - 9 ) 
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(4.1-10) 



The float time constant Tf is computed by considering only the gimbal 
moment of inertia about the output axis and output axis damping: 


T i 



(4.1-11) 


Since only the linear parts of the loop are changed as a result 
of considering compliance, many of the phenomena seen for simple 
gyro dynamics are again observed. As in Section 4. 1. 1, the effect of the 
sample and hold operation is neglected until later in this discussion. Due 
to the absence of any phase shift contribution from the nonlinearity, the 
loop still limit cycles at a frequency determined only by the phase char- 
acteristics of the linear part of the open loop. As before, the limit cycle 
"linearizes" the nonlinearity with respect to a signal circulating in the 
loop. Unfortunately, the closed loop signal transfer characteristic can 
no longer be approximated by that of a second order system. 

The limit cycle analysis of the loop under discussion is more 
difficult than when simple gyro dynamics are considered. Here, the linear 
open loop transfer function H(s) has seven poles, instead of three. 

Ktg / \ 

H(s) = sr^Tl\M^ t ) K sg (4.1-12) 
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where C4, C3, C2, Cj and C 0 are defined above. It is observed that the 
phase condition for a limit cycle, £H(jco) = -180°, occurs at several fre- 
quencies. The limit cycle frequencies can be computed numerically as 
the positive real root(s) of a polynomial which is derived from this con- 
straint. 

y 10 w 10 + y g w 8 + y 0 w 6 + y 4 w 4 + y^co 2 + y Q = 0 (4. 1-14) 
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In general, there can be as many as five positive, real solutions to Eq. 

(4. 1-14), each corresponding to a different "crossover frequency” where 
the linear, open loop phase lag is 180°. Stability analyses are applied to 
discard certain solutions of Eq. (4. 1-14). In general, if the five solutions 
are ordered by ascending frequency (oo^, u^, ° * *, o^g), the stable limit 
cycles correspond to the frequencies with odd numbered subscripts (co^, 

For each stable limit cycle there is a different limit cycle amp- 
plitude, corresponding to a different gain, N B , of the nonlinearity to signal 
inputs. Because the closed loop frequency response to an input depends on 
N b there can be as many as three different sets of transfer characteristics, 
each resulting from a different limit cycle. Figure 4. 1-5 shows two closed- 
loop gain responses obtained for a system with two stable limit cycle modes 
(coj % = 727 Hz, 868 Hz). Although the limit cycle frequencies differ signifi- 
cantly, the two responses are quite similar. Appendix G discusses the fact 
that exciting the gyre loop with a particular input at certain frequencies can 
cause the limit cycle to change , with a resultant shift to a different frequency 
response. The system whose closed loop is shown in Fig. 4. 1-5 can be 
shifted from the 727 Hz limit cycle mode to the 868 Hz mode by excitation 
with a 0. 1 rad/sec amplitude angular rate input at 810 Hz or above. How- 
ever, it cannot be shifted back to the 727 Hz mode in a corresponding manner. 

Two effects - finite sampling frequency and limit cycle quench- 
ing - remain to be discussed in connection with binary loops containing 

t 

rotor- to-gimbal compliance. It was shown in Section 4. 1. 1, that the 
sample and hold operation introduces a lag which can be approximated by 
the transfer function 1/(1 + jvT/2) where T is the sampler period. This 
causes a modification of the limit cycle frequency and amplitude, changing 
the nonlinearity gain to the signal and altering the closed loop response. 
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Closed Loop Gain 



Figure 4.1-5 Potential Frequency Responses of 

System with Two Stable Limit Cycles 
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The effect of sampling on a gyro with significant compliance between rotor 
and gimbal is illustrated in Fig. 4. 1-6. Two analytically determined fre- 
quency response curves are shown; one for a high sampling frequency 
(T = 0) and the other for a 10 kHz rate. Limit cycle quenching and its 
effect on the performance of the gyro loop discussed in this section are 
treated in Appendix D. 


4. 2 SIMULATION OF THE LIMIT CYCLING BINARY GYRO 

The analyses of Section 4. 1 provide interesting insights into 
the closed loop behavior of the binary torqued gyro. The accuracy of 
these predictions depends on how well the describing function linearization 
characterizes the true nonlinearity response to a signal in the presence 
of a limit cycle. As a check on this approach, typical gyro loops were 
analyzed by the methods discussed in Section 4. 1 and were then simulated 
on an analog computer to verify the gyro transfer characteristics predicted 
analytically. 


4.2.1 Simple Gyro Dynamics 


A binary gyro loop was examined analytically and by analog 
simulation. The dynamic characteristics of the linear portion of the gyro 
were identical with those presented in Fig. 4. 1-1. Details of the simula- 
tion are given in Appendix F. The important gyro parameters for this 
comparison are: 


H = 2 x ufi dyne-cm-sec 

max = lrad/sec 
Tf = 0. 0025 sec 

T tg = 0.0001 sec 
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Frequency (Hz) 


Figure 4.1-6 


Effect of Sampling Rate on Gyro 
Frequency Response: Significant 
Rotor Bearing Compliance 





Figure 4.2-1 displays simulation records of the binary gyro 
float angle and torquer output in the absence of an input angular rate. A 
very high sampling rate (1 mHz) was used. The limit cycle observed was 
in close agreement with that predicted for a continuous loop. 



(b) Float Angle, « 0 


Figure 4.2-1 Simulation Records for Limit-Cycling 

Binary Gyro 
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An oscillatory angular rate with a peak amplitude of 0. 1 rad/sec 
was applied to the gyro loop. Figure 4.2-2 and 4.2-3 display simulation 
records of gyro input, float angle and torquer output for two values of input 
frequency. Again, an essentially continuous loop was considered. It can 
be seen that for a low frequency input the gyro float angle follows the input 
very closely, despite the simultaneous presence of a limit cycle oscillation. 
At a higher input frequency, however, it is more difficult to identify signal 
and limit cycle effects separately. The closed loop frequency response 
predicted analytically for this system is seen in Fig. 4.2-4, together with 
points obtained by simulation. There is a good agreement between the two 
techniques. 

Figure D-2 of Appendix D indicates that, for an input amplitude 
of 0. 1 rad/sec, the limit cycle seen in this system is quenched over a 
range of frequencies above 200 Hz. Accordingly, the analytically derived 
response has been corrected over this frequency range to show a flat- 
topped peak characteristic cf the binary system when the limit cycle is 
quenched. In this region the system output is a square wave at the input 
frequency with amplitude equal to the nonlinearity drive level. This level 
is equivalent to 1 rad/sec input rate and the fundamental component of the 
output has an amplitude of 4/n = 1. 273 rad/sec. Since the input amplitude 
is 0. 1 rad/sec, gain of the closed loop is a constant value of 12.73. More- 
over, all information about the input signal amplitude (except that it is 
large enough to extinguish the limit cycle) is lost, and the gyro loop ceases 
to function as a linear system. Limit cycle quenching can represent a 
definite restriction on the range of the amplitudes and frequencies that can 
be accurately detected by a binary gyro. However, quenching is not com- 
mon at frequencies well below that of the limit cycle. 
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(b) Float Angle, a D 



(c) Torque r Output 

Figure 4.2-2 Binary Gyro Simulation: 20 Hz Input 
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(a) Input Angular Rate, <4 



(b) Float Angle, a Q 



(c) Torquer Output 

Figure 4.2-3 Binary Gyro Simulation: 165 Hz Input 


Sampling - Simulations were also conducted for the binary limit 
cycling gyro with significant sampling lags. Figure 4.2-5 shows the fre- 
quency response for a loop with 10^ samples per second. It differs notice- 
ably from that shown in Fig. 4. 2-4. The solid line indicates the response 
computed using describing function theory and simulation points obtained 
are also shown. Good agreement between analysis and simulation is indi- 
cated. 
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Gain by Simulation 
Phase Shift by Simulation 
Gain Analysis 
Phase Shift Analysis 

Infinite Sampling Rate 



Frequency (Hz) 

Figure 4. 2-4 

Frequency Response of a Limit Cycling 
Binary Gyro: Comparison of Analytic 
and Simulation Results 





Frequency (Hz) 


Figure 4.2-5 


Frequency Response of a Limit Cycling 
Binary Gyro: Comparison of Analytic 
and Simulation Results 





4.2.2 Gyro Dynamics with Rotor Bearing Compliance 

Section 4. 1.2 discusses the changes made to the gyro model 
when rotor- to- gimbal compliance is assumed to be significant. For the 
purpose of comparing analytic and simulation results, the following set 
of gyro parameters was employed: 

Wg = 4000 rad/sec 
v n = 2512 rad/sec 

p =0.02 

The frequency response predicted analytically for this system is shown in 
Fig. 4.2-6, along with points obtained by simulation. Good agreement is 
observed between the two methods except near the frequencies at which 
limit cycle quenching takes place. A small hysterisis region, not pre- 
dicted by the analysis, is observed in the simulation for input frequencies 
near the first (low frequency) quenching region. The other quenching 
region, occurring at frequencies near that of the limit cycle, is not accur- 
ately described by the analysis. This error in determining quenching is 
thought to be caused by beat frequency terms not considered in the describ- 
ing function approach used. It can be seen from the figure that the analysis 
provides adequate estimates of the general frequency range for quenching 
and gives an excellent prediction of the low frequency portion of the gyro 
loop response. Comparison of Figs. 4. 2-4 and 4. 2-6 indicates that, for 
the parameters chosen, rotor bearing compliance reduces the useful 
bandwidth of the gyro. Since the low frequency peak in Fig. 4.2-6 is 
caused by the compliance, greater bearing stiffness would help to alle- 
viate this effect. 
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Input Frequency (Hz) 


Figure 4. 2-6 


Frequency Response of a Binary Gyro with 
Rotor Bearing Compliance: Comparison 
of Analytic and Simulation Results 






4.3 COMPENSATION FOR THE LIMIT CYCLING BINARY GYRO 

Introduction - The binary gyro loop limit cycles whether or not 
a rate signal is present at the input (unless the signal is of high enough 
amplitude or suitably shaped to quench the limit cycle). As a consequence 
of the limit cycling property, the closed loop frequency response to an 
input signal is largely determined by gyro and torquer time constants and 
not by easily manipulated electrical gains (see Section 4. 1. 1). Specifically, 
if the loop is linearized using describing function techniques, a response 
characteristic approximately that of a second order system can result with 
natural frequency and damping ratio given by Eqs. (4. 1-1) and (4. 1-2). 
These frequency response characteristics cannot be changed without making 
hardware modifications to change and \ g. In general, Tjg is made as 
small as possible in order to permit high sampling rates. Also the float 
time constant is related to gyro noise effects through the damping coeffi- 
cient, C. Raising the damping reduces gyro resolution and increases Hie 
effect of electrical noise in the pickoff. In addition, no reasonable choice 
of the time constants can make the damping ratio greater than 1/2. Thus, 
the frequency response is forced to have a peak near the natural frequency. 
Clearly, some frequency sensitive compensation may be desired to remove 
this peak and, more generally, to permit control of the frequency response 
characteristics without changing Tj and r^g. 

Analysis - Following the approach used to obtain Eqs. (4. 1-1) 
and (4. 1-2), we represent the effective gain of the binary nonlinearity to 
a sinusoidal signal by the Two Sinusoid Input Describing Function (TSIDF), 
written Ng(A,B). This describing function is the ratio of the amplitude of 
the signal- frequency harmonic at the nonlinearity output to the amplitude 
B of the signal at the nonlinearity input when a sinusoid (limit cycle) of 
amplitude A is also present at the input. If B is considered small com- 
pared to A, a simplified form for Ng(A,B) results: 
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N B (A,B) * N b (A,0) B « A 

(4.3-1) 

a 2D 
irA 

where D is the drive level of the nonlinearity. Appendix C examines in 
greater detail the conditions under which this approximation holds. 

In attempting to compensate the limit cycling binary loop, we 
must consider the effect of the limit cycle on the total open loop gain and 
phase performance as seen by the signal sinusoid. The binary nonlinearity 
TSIDF gain to the limit cycle, written N A (A, B), is given by 

N A (A, B ) = |j; B«A (4.3-2) 

The complex open loop transfer function of the linear part of the loop is 
written H(jw). The definition of the limit cycle requires that, at the limit 
cycle frequency, city, 


H(joty) N A (A,B) = 1 A 180° (4.3-3) 

The limit cycle amplitude A and frequency adapt to maintain this equal- 
ity regardless of the compensation introduced. From Eqs. (4.3-1) and 
(4.3-2), it can be seen that 

N A (A,B) = 2N B (A,B); B « A (4.3-4) 

Substituting for N A in Eq. (4.3-3) yields 


H(jw A )N B (A,B) = j A1S0° (4.3-5) 
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Here, H(jWje) Ng(A,B) is the expression for the complex open loop gain as 
seen by the signal sinusoid. We conclude that, due to the adaptive proper- 
ties of the limit cycle amplitude and frequency, the open loop gain and 
phase shift as seen by the signal sinusoid will always be 1/2 /- 180° at coj. 
The linear analysis technique used to investigate in- loop compensation of 
binary pulse torqued gyros makes use of the rectangular gain-phase plot 
(Nichols chart). The open loop gain and phase characteristics are plotted 
as a function of frequency. On the same plot, lines of constant M (closed 
loop amplitude ratio) are constructed. The closed loop amplitude ratio 
for any sinusoidal input is the value of M intersected by the open loop 
gain-phase curve at the same frequency. The M = 1 line represents the 
desired unity closed loop gain. Figure 4.3-1 shows the M = 1 line and the 
open loop gain-phase curve of the uncompensated gyro loop with the follow- 
ing parameters: 


= 10 “^ sec 
r f • =' 2.5 x 10"3 sec 

Note that the gain-phase curve cuts through high closed loop gain lines at 
frequencies lower than 004 , producing a peaked frequency response. 

In Section 4. 1 torquer lag was neglected to allow approximation 
of the uncompensated binary gyro loop by a second order linear system. 
This representation of the gyro is adequate for computing closed loop gain 
and phase shift characteristics at low frequencies. However, the torquer 
begins to contribute open loop phase shift at frequencies below the limit 
cycle, causing a significant difference between the open loop gain-phase 
plots of the approximate (second order) and more complete gyro models. 
Consequently, the peak closed loop gain calculated using a second order 
approximation to the gyro loop is always less than the true value. 
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Figure 4. 3-1 Open Loop Gain- Phase Plot for 

Uncompensated Binary Gyro 


In order to obtain unity closed loop gain at all frequencies below 
the limit cycle, we must make gain and phase changes to the open loop 
which cause the compensated gain- phase curve to coincide with the M = 1 
line. As discussed earlier, the complex open loop gain must always be 
1/2 7-180° at the limit cycle frequency, coj; the gain- phase curve always 
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crosses the -180° phase shift line with a gain of 1/2 at whether or not 
compensation is used. It can be seen from Fig. 4.3-1 that the M = 1 line 
also intersects the -180° line at 1/2. Four sets of linear dynamic compen- 
sation schemes are discussed below. In each case the compensation is 
assumed to be placed in the forward signal path as illustrated in Fig. 4. 3-2. 
The closed loop gains for gyros employing various compensation schemes 
are shown in Fig. 4.3-3. 



Figure 4. 3-2 Position of Binary Gyro 

Loop Compensation 


From Fig. 4. 3- 1 it can be seen that the linear compensation 
necessary to provide an open loop gain- phase characteristic which coin- 
cides with the unity closed loop gain curve must provide rapidly increasing 
phase lag with little change in amplitude ratio at frequencies immediately 
below cojj. A resonant second order system provides rapid changes in phase 
shift but gives a sharp amplitude ratio peak. By cascading such a system 
with a first order lag to reduce the gain peak, it was hoped that the desired 
gyro open loop characteristics could be achieved. For the gyro parameters 
postulated, this form of compensation produces both amplification and at- 
tenuation in the critical range of frequencies, <o< cty. Departure from 
unity closed loop gain is not desirable at any frequency. Large amplification 
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in particular is associated with limit cycle quenching. In view of the fact 
that this compensation scheme cannot produce a perfectly flat response, 
the best course of action is to choose a set of compensation parameters 
which will minimize the maximum deviation of any part of the closed loop 
response from unity gain. The values used to get the response shown are: 


resonant natural frequency = 159 Hz 

resonant damping ratio = 0. 03 

lag time constant = 3.3 x 10" 4 sec 

The resulting closed loop response has a maximum gain of approximately 
1.29 at 48 Hz and a minimum gain of about 0. 80 at 100 Hz with a limit 
cycle frequency of 126 Hz. (The uncompensated gyro gain has a peak of 
about 7. 5 at 230 Hz and Wj = 318 Hz. ) The decrease in limit cycle frequency 
as a result of compensation is unavoidable. Since this form of compensa- 
tion adds open loop phase lag, is lowered and the unstable bandwidth is 
always decreased when the binary limit cycling loop is compensated in this way. 
Also, because of the nature of binary torqued gyro loops, the limit cycle 
amplitude increases. 

Another attempt to reduce variations in the closed loop gain 
uses a lead- lag- lag- lead network to generate the abrupt open loop phase 
shift required. The network transfer function is 


(saT c + 1) (s — + 1) 


(ST- + 


W 


a > 1 


The parameter a is chosen to be 20. Selection of T c = 1/toj causes the 
compensation to produce no phase shift at Wj. Thus, the limit cycle fre- 
quency is unchanged. 
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Closed Loop ( lain 



Figure 4. 3-3 Closed Loop Gain for Limit Cycling 

Binary Gyro: Effect of Loop 
Compensation 


If it is desired to increase the useful bandwidth of the loop with- 
out attempting to significantly alter the peak gain, lead- lag compensation 
is effective. A compensation network with the following transfer function 
was used to get the response illustrated in Fig. 4.3-3. 

sar c + 1 a = 20 

sr c + 1 ’ r c = 0. 00005 sec 
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Integral-Bypass Compensation - The motivation for using inte- 
gral bypass compensation is to reduce the average low frequency float 
angle (taken over one or more limit cycle periods). Specifically, it is 
desired to use compensation at the nonlinearity input which approximates 
a pure integrator for low frequencies, thus giving the open loop a double 
integration characteristic for constant inputs. When the additional inte- 
gration is inserted in the loop, the average float angle is zero for a con- 
stant input. This compensation also reduces the float angles caused by 
low frequency inputs. 



o 


Figure 4. 3-4 Block Diagram for Integral- 

Bypass Compensation 


Integral-bypass compensation (see Fig. 4. 3-4) has the transfer 

function 


s + K 
s 

At low frequencies the transfer function is approximately K/s, while at 
high frequencies the compensation has unity gain, no phase shift and does 
not influence the loop response. The effect of integral bypass compensa- 
tion is seen in Fig. 4. 3-5 as a bending of the open loop response towards 
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Figure 4.3-5 Nichols Chart Showing the Effect of 
Integral Bypass Compensation on 
Open Loop Response 
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the -180° phase line at low frequencies. Increasing K broadens the 
frequency range over which float motion is suppressed by enlarging the 
frequency band where the compensation behaves like an integrator. 

Note that a value of K as high as 10 has little effect on the limit cycle 
frequency or the closed loop response near a > 4 . For values of K s 10, 
integral bypass compensation reduces the average float hangoff from low 
frequency inputs and leaves the closed loop response and limit cycle es- 
sentially unchanged. 

Conclusions - The value of compensation in the limit cycling 
binary gyro loop has been studied analytically. Two forms of compensa- 
tion which reduced the peak closed loop gain gave poor loop behavior at 
low frequencies. A third form of compensation, introducing lead into the 
gyro loop, increased the bandwidth of the gyro loop. Integral compensa- 
tion was investigated as a means for stiffening the response of the gyro 
float to low frequency inputs. The investigation conducted was exploratory 
in nature and no consideration was given to the practical aspects of im- 
plementing the compensation schemes proposed. 


4.4 LINEAR ANALYSIS OF THE DITHERED BINARY GYRO 

The addition of a dither signal to the nonlinearity input of a 
binary torqued gyro can be used to create an approximately linear relation 
between the gimbal angle, a 0 , and the gyro output (see Ref. 1). Also, if 
attenuation of the dither oscillation by float dynamics is large enough, the 
dither nonlinearity can be represented as a limiter as shown in Fig. 4.4-1. 
Linear feedback analysis techniques can be used to investigate the transfer 
characteristics of the loop illustrated over a wide range of amplitudes and 
frequencies. However, the limiting nature of the new representation for 
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the dithered binary nonlinearity raises the possibility of a no- input limit 
cycle in the rebalance loop. Careful consideration reveals that the poten - 
tial limit cycle in the dithered gyro is essentially the same unforced 
oscillation that occurs in any simple binary loop . Using the block diagram 
of Fig. 4.4-1, it can be seen that any limit cycle must have a frequency 
which provides it radians of phase lag in the linear part of the loop; the 
nonlinearity shown provides no phase change. The sinusoidal describing 
function for the limiter, N^, is illustrated in Fig. 4.4-2. Peak gain 
occurs when the input sinusoid has an amplitude smaller than the dither 
magnitude, A<j. In that case the describing function is a constant, D/A^, 
where D is the torquer current level. For a limit cycle to exist, the fol- 
lowing relation must be satisfied: 

N A |H(jWi)| = 1 (4.4-1) 

For the gyro model illustrated the open loop gain at the limit cycle fre- 
quency is; 


ft'*)! = (4.4-2) 

C r f + r tg 

Consequently, to prevent a limit cycle the dither amplitude must be chosen 
so that 


KggKjgD 


T f T tg 
T f + T tg 


(4.4-3) 


If the dither frequency attenuation by the float is sufficient (validating the 
representation of the nonlinearity shown in Fig. 4.4-1) and Eq. (4.4-3) is 
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satisfied, the application of dither at the input to the binary nonlinearity 
will prevent the occurrence of a limit cycle. Analysis of the dithered gyro 
loop for a wide range of inputs can then proceed using techniques developed 
for linear feedback systems . 


Float Dynamics 



Figure 4.4-1 Dithered Binary Gyro 


"a 



Figure 4.4-2 Describing Function for the Limiter 
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Example - Consider a single-degree-of-freedom gyro with the 
following parameters: 

Iq 0 = 250 gm-cm^ 

Kjg = 1200 dyne-cm/mA 
H = 2 x 10® gm-cmVsec 
C = 1 x 10® dyne-cm-sec 
K sg = 20 mv/mrad 

“Vax" DK tg/ H = 1 rad / sec 

Using Eq. (4.4-3), the minimum dither amplitude to preclude a limit 
cycle is 3.84 mV and the slope of the linear portion of the nonlinearity, 
Kj = D/A d , is 43.2 mA/mV. The natural frequency and damping ratio of 
the linearized gyro loop (ignoring torque r dynamics) are 


= •y/ Kss ^ t g* ^ = 325 Hz 


C = 


1 

2T fWn 


0.098 


(4.4-4) 


Figure 4.4-3 illustrates the resultant frequency response for the linearized 
loop. For purposes of comparison the frequency response of a limit cycling 
binary gyro with the same parameters is shown. In interpreting the curve 
shown it should be remembered that for the dithered binary gyro no oscilla- 
tion faster than half the dither frequency can be accurately detected at the 
gyro output. This is a consequence of averaging over each dither cycle. 

As the dither amplitude is increased above the minimum value required 
to prevent a limit cycle, Kj decreases. From Eq. (4.4-4) it can be seen 
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that the natural frequency of the linearized gyro decreases and its damping 
ratio increases when is made larger,, The amplitude ratio frequency 
response for the dithered gyro is also shown in Fig. 4.4-3 for the case 
where dither amplitude is twenty-five times the minimum value computed 
above. Note the difference from that achieved by the limit cycling gyro. 



Figure 4.4-3 Frequency Response Characteristics 

for Binary Rebalance Loop 


When the dithered binary gyro is represented as illustrated in 
Fig. 4.4-1 the frequency response characteristics of the second order 
linearized rebalance loop can be controlled by varying dither amplitude. 
However, damping ratio and natural frequency cannot be controlled 
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independently. Any increase in damping ratio (increase in dither ampli- 
tude) necessarily implies a decrease in natural frequency. In the example, 
the entire range of desirable loop transfer characteristics (co n and C ) can 
be generated using dither signals large enough to prevent a limit cycle. 
Considerably more design flexibility is available with the dithered binary 
gyro. In addition, the interdependence of peak closed loop gain and band- 
width can be eliminated through insertion of linear compensation in the 
loop. The whole range of linear systems techniques is available for ana- 
lyzing the dithered loop. 

The use of a dither signal with a ternary nonlinearity offers no 
advantage relative to the technique discussed above. 
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5 . 


TERNARY PULSE REBALANCED GYRO 


SUMMARY 


Research Performed - The behavior of ternary pulse 
torqued gyros is predicted analytically. A major por- 
tion of the effort is devoted to determining the closed 
loop gain characteristics of ternary gyros and verify- 
ing the analysis by simulation. The influence of rotor 
bearing compliance is also analyzed and the predicted 
gyro response is checked by simulation. The effects 
on the ternary gyro frequency response of different 
gyro input amplitudes and random noise at the logic 
inputs are also treated. 

Results - It is demonstrated that the closed loop be- 
havior of ternary pulse rebalanced gyros can be accur- 
ately predicted. The gain and phase characteristics 
observed during simulation are in close agreement with 
those computed analytically. 

It is shown that the ternary pulse rebalanced gyro can 
exhibit a jump resonance in its closed loop response to 
sinusoidal inputs. This history- dependent gain and 
frequency shift is observed in both the analysis and the 
simulation. In addition, analysis indicates that the gyro 
has a transfer characteristic which can depend greatly 
on the amplitude of the input signal; the ternary gyro 
loop transmission is nonlinear. The latter property 
was not thoroughly checked by simulation but enough 
confidence has been gained in the analytic technique that 
the conclusion is considered to be correct. 

The effect of rotor bearing compliance on the ternary 
gyro loop response is similar to that found for the bi- 
nary gyro. Both analysis and simulation indicate that 
bearing compliance can reduce the useful bandwidth of 
the gyro. Good agreement between the two approaches 
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is observed, verifying the utility of the analytic techni- 
que. The addition of Gaussian noise to the signal rep- 
resenting float angle is analytically shown to create a 
more nearly linear transfer characteristic in the ternary 
gyro. 


Conclusions Related to Gyro Design - It is demonstrated 
that the characteristics of ternary pulse rebalanced 
gyros can be predicted analytically. However, consid- 
erably more effort is required than for the binary gyro. 
No simple set of equations permits quick calculation of 
ternary gyro response. There is no limitation on the 
complexity of the linear gyro loop dynamics which can 
be treated by this approach. 

Non-unique closed loop gains are predicted and observed 
for the ternary gyro, indicating that this instrument 
can exhibit an ambiguous response for inputs in a cer- 
tain frequency range. In addition, analysis indicates 
that in a large band of input frequencies the ternary 
gyro loop gain depends on the input amplitude. These 
two observations suggest that a ternary pulse rebalanced 
gyro should not be used where high frequency angular 
vibrations must be measured accurately. As with the 
binary gyro, rotor bearing compliance can reduce the 
bandwidth of the ternary gyro. 

Analysis indicates that the presence of random noise 
in the signal representing float angle can cause the 
ternary gyro to behave in a more linear fashion. How- 
ever, the technique employed in the analysis does not 
permit phase shift characteristics to be specified. 
Further investigation is required before a random 
noise (dither) is intentionally added in the signal elec- 
tronics of a ternary gyro to give a better response. 
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Introduction - Describing function analysis was shown to be of 
considerable value for treating the binary gyro in Section 4. Similar 
techniques are applied below to the ternary gyro of Fig. 5. 1-1. Though 
the approach is basically the same, the effort involved is greater. Again, 
in order to simplify the problem the sampling nature of the pulsed gyro, 
loop is ignored. 


u>. 

1 



A 

W i 


6= dead -zone threshold in mV 
D = drive level in mA 

Figure 5. 1-1 Gyro Pulse Rebalance Loop 

with a Ternary Nonlinearity 


5. 1 ANALYSIS OF THE FREQUENCY RESPONSE 
CHARACTERISTICS OF THE TERNARY GYRO 

As in Section 4, we desire to make approximations which will 
permit application of linear analysis techniques for determining the closed 
loop response to sinusoidal inputs. With the binary gyro, the presence of 
a limit cycle allowed us to approximate the describing function gain to a 
signal sinusoid by a quantity which depends only on the limit cycle amplitude 
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The closed loop frequency response to the signal was then determined by re- 
placing the nonlinearity with the appropriate gain and applying conventional 
linear analysis techniques. Unfortunately, the ternary gyro does not exhibit 
a clearly defined limit cycle and the describing function gain of the three 
level nonlinearity must be treated as a function of the amplitude of the signal 
entering it. Also, because the float provides frequency dependent attenua- 
tion, the amplitude of a sinusoid entering the nonlinearity is a function of 
both the magnitude and frequency of the gyro loop input. Consequently, a 
different size sinusoid appears at the input of the nonlinearity for each 
frequency and amplitude combination applied at the loop input. For 
each input combination explored in the course of a frequency response ana- 
lysis, a different nonlinearity equivalent gain must be calculated and a sep- 
arate determination of system output amplitude made- -all to obtain one 
value of closed loop gain. To simplify matters in the example which follows, 
the input amplitude is held constant and only the frequency varied to obtain 
a plot of closed loop gain. However, it should be emphasized that a different 
curve will result from each value of input amplitude assumed. 

Linearization of the Ternary Nonlinearity - The gain approxima- 
tion used to describe the ternary nonlinearity is the Sinusoidal Input Des- 
cribing Function (SIDF), written N^(A) and defined as the ratio between the 
amplitude of the fundamental sinusoid at the nonlinearity output and the 
amplitude A of the pure sinusoidal input to the nonlinearity. N^(A) is 
evaluated from: 

n a(A) = Is Vl - (j) (5.1-1) 

It should be noted that the SIDF gain relates only to the first harmonic of 
the waveform at the nonlinearity output. All higher harmonics are ignored. 
They are assumed to be attenuated when passed through torquer and float 


80 



dynamics. The assumption that the higher harmonics are removed by low 
pass elements in the loop is generally correct when treating that portion of 
the output signal reaching the nonlinearity input. Consequently, the closed 
loop analysis is not invalidated by ignoring higher harmonics at the output. 
However, high frequency signals do exist at the gyro output and will enter 
the direction cosine calculations. The frequency response curves obtained 
in this section account only for the first harmonic of the gyro loop output. 

The signal entering the nonlinearity is essentially the sum of the output funda 
mental and the loop input signal, both filtered by float dynamics. The 
sampling nature of the nonlinearity is ignored in the analysis which follows. 

Analysis Example: Simple Gyro Dynamics - The gyro loop is 
divided into two parts (see Fig. 5. 1-2) each of which will be treated sep- 
arately to obtain certain characteristics. Then the two similar sets of re- 
lations are compared to obtain a simultaneous solution for the system. The 
technique is analogous to load line solutions in electrical engineering. Linea 
analysis is applied to the portion of the system shown in Fig. 5. l-2a. If 
cq is a sinusoid of known magnitude, an amplitude relation can be obtained 
between two sinusoids of the same frequency as the input, one at point a 
and the other at point b. This ratio results from linear theory and does not 
depend on the specification of an open loop gain relation between the signals 
e a and e^ (see Appendix E). When plotted for a particular gyro, the relation 
defines an ellipse for each different amplitude and frequency of the input 
sinusoid. If the input magnitude is assumed constant, a set of ellipses, 
illustrated in Fig. 5.1-3, results. Each ellipse represents a different input 
frequency. 


The ternary nonlinear element is treated separately (see Fig. 

5. l-2b). The describing function, Eq. (5. 1-1), is used to obtain a relation 
similar to that computed from linear system theory. This characteristic, 
illustrated in Fig. 5. 1-4, depends on the quantities D and 6. When the plots 
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representing the two parts of the loop are superimposed in Fig. 5. 1-5, 
the intersections constitute a set of operating points for the system. 

Since each ellipse represents a different loop input frequency, each inter- 
section with an ellipse represents a loop operating point at that frequency. 
Every intersection point also supplies a particular value for the nonlinearity 
output. Since the output of the nonlinearity is also the loop output, the 
absolute value of the loop transfer function is simply e^ divided by the input 
amplitude. The ratios obtained for the ellipse of a particular frequency, 
v, represent possible closed loop system gains to a sinusoid at that fre- 
quency. However, not all of the intersections represent stable solutions. 



Figure 5. 1-4 Output vs Input Computed from 

the SIDF Gain for the Nonlinearity 


By determining a sufficient number of intersection points, it is 
possible to construct a closed loop amplitude ratio response curve for the 
ternary gyro. Figure 5. 1-6 illustrates such a curve drawn from the inter- 
section points determined in Fig. 5. 1-5. The peaked portion of the closed 
loop response curve exhibits a ’’jump resonance” characteristic similar to 
that given by a softening spring. The effect of deadzone on the closed loop 
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gain is seen as an abrupt drop to zero at the frequency v c . If there was 
no deadzone , the response gain would fall off gradually with increasing 
frequency. However, the deadzone causes the relay to block all signals 
whose amplitude at the nonlinearity input is less than 6. As frequency in- 
creases, the float gain decreases until, at frequency v c , input amplitude 
at the nonlinearity input is less than 6. Then the output, and consequently 
the loop gain, drops abruptly to zero. Deadzone also causes the overhang- 
ing nature of the response curve. At high input frequencies, attenuation 
by float dynamics prevents the input to the nonlinearity from exceeding 6. 

As frequency is decreased the float angle eventually becomes larger than 
the deadzone, providing a loop output and generating signals at all points in 
the loop. This occurs at a frequency different from v c because the gyro 
is essentially an open loop system until an output appears at the nonlinearity. 



Figure 5.1-5 Superposition of Figs. 5.1-3 and 

5. 1-4 Showing Intersection Points 
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Closed 

Loop 

Gain 



Figure 5„ 1-6 Ternary Gyro Closed Loop Gain 


Effect of Input Amplitude - Because the nonlinearity it contains 
has an amplitude dependent output, the ternary gyro loop exhibits amplitude 
dependent characteristics. When plotted for different magnitudes of the 
input sinusoid, the frequency response curves have various shapes (see 
Fig. 5. 1-7). If the input amplitude is small enough (0.01 rad/sec in the 
figure) the closed loop response has no peak. As the size of the input 
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increases the jump resonance characteristics and dual gain cutoffs dis- 
cussed above become more pronounced. A look at Fig. 5. 1-5 provides 
an explanation. The ellipsoidal curves obtained from linear analysis are 
changing size and rotating. As frequency increases, the length of the 
major axis first increases, then declines. However, rotation of the major 
axis is monotonic. For the smallest input amplitude plotted in Fig. 5.1-7 
no appreciable rotation of the ellipses takes place before the loop output 
drops to zero. When the gyro input is raised to 0.03 rad/sec the ellipses 
rotate and the intersection points occur at higher values of e^, giving a 
peaked loop response. When the largest input is analyzed, the ellipses 
not only rotate but also become long enough to provide several intersection 
points at some frequencies. A jump resonance characteristic results. 

The fact that ternary gyros have a gain which depends on the input amplitude 
is a serious disadvantage for this type of instrument. If it is necessary 
to accurately measure high frequency angular vibrations, the ternary pulse 
rebalanced gyro does not appear satisfactory. 

Gyro Dynamics when Rotor Bearing Compliance is Considered - 
The loop with rotor- to- gimbal compliance is treated in the same manner 
as that with simple dynamics discussed above. Ellipses are generated to 
describe the linear part of the loop using the equations in Section 4. 1. 2. 

The SIDF gain for the nonlinearity remains unchanged. Points of inter- 
section between the ellipses and the nonlinearity input- output characteris- 
tics define the loop output amplitude as a function of frequency. From 
this information, the closed loop gain to a sinusoid is inferred. Consider- 
ation of rotor bearing compliance does not change the basic procedures 
followed in obtaining the closed loop frequency response from describing 
function analysis. 
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5.2 SIMULATION OF THE TERNARY GYRO 


The analytic techniques demonstrated in Section 5. 1 provide 
interesting insights into the closed loop behavior of the ternary-torqued 
gyro loop. However, the accuracy of these predictions depends on how 
faithfully the describing function approximation characterizes true non- 
linearity performance. As a check on the validity of this approach, char- 
acteristics of certain ternary loops were computed and the results verified 
using analog computer simulation. 


5.2.1 Simple Gyro Dynamics 


A set of parameters is chosen to 
gyro loop with second order gyro dynamics 


represent a typical ternary 
(see Fig. 5. 1-rl): 


“W ■ 1 rad/sec T tg 

Tj = 0. 00125 sec 6 

H/C = 1 


0. 00005 sec 

4. 17 x 10"^ rad about the 
input axis 


With the exception of the nonlinearity deadzone 6, these are the same 
values as those shown in Fig. 5.1-7. A different value of deadzone was 
used for this comparison to ensure that no limit cycle took place in the 
simulated gyro loop. Oscillatory angular rate inputs with a peak amplitude 
of 0. 1 rad/sec were applied to the simulated gyro. Figures 5.2-1 and 
5.2-2 show records of the input, float angle and torquer output when two 
different input oscillation frequencies were simulated. The sampling rate 
in both cases illustrated is 1 mHz providing essentially continuous ternary- 
loops. The closed loop frequency response of this gyro is distinguished by 
a peak at 400 Hz (see Fig. 5.2-3). Discrete points obtained via simulation 
are indicated while the curve predicted by describing function analysis 


88 




(b) Float Angle, a 0 



(c) Torquer Output 

Figure 5.2-1 Ternary Gyro Simulation Records: 

20 Hz Input 
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(a) Input Angular Rate, 



(b) Float Angle, a 0 



(c) Torquer Output 


Figure 5.2-2 Ternary Gyro Simulation Records: 

150 Hz Input 


0 


appears as a solid or dashed line. Note that the analytic response contains 
unstable segments in the jump resonance region which cannot be observed 
for any length of time by simulation. They represent unstable modes which 
do not occur in the actual instrument either. Agreement between the analy- 
tic and simulation results is quite good on the stable portions. The response 
is similar in shape to the curve seen in Fig. 5. 1-7. 



Figure 5. 2-3 Frequency Response of a Ternary Gyro: 

Comparison of Analytic and Simulation 
Results 
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Figure 5. 2-4 presents simulation records illustrating the behavior 
of the ternary gyro float angle at a jump resonance point. The input fre- 
quency is being raised so slowly that the change cannot be detected on the 
record, but the gyro float angle behavior shows a clear discontinuity. The 
gain of the gyro loop was observed to increase suddenly at the point where 
a 0 becomes larger. 
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0. 1 mrad 


(b) Float Angle, <x 0 


Figure 5. 2-4 Ternary Gyro Simulation Record 

Illustrating Jump Resonance 




5.2.2 Gyro Dynamics with Rotor- to- Gimbal Compliance 


The same gyro discussed in Section 5. 2. 1 is treated here for 
the case where compliance between the rotor and gimbal are significant. 

The parameters related to compliance are: 

tog = 4000 rad/sec 
= 2512 rad/sec 

p =0.02 

The float time constant, T^, is not the same quantity when the rotor and 
gimbal are not rigidly connected. The value = 0. 001 sec was used. The 
frequency response curve for this system has a peak centered near 180 Hz 
as shown in Fig. 5.2-5. The solid line represents the analytic solution and 
simulation results appear as discrete points. As before, certain points in- 
dicated analytically are unstable and cannot be observed. The simulation 
points substantiate quite well the response obtained by describing function 
analysis, however. 

Note that the response peak resulting from consideration of com- 
pliance dynamics occurs at a noticeably lower frequency than the peak asso- 
ciated with a second-order description of the gyro. Apparently, the closed 
loop poles are radically altered by including compliance dynamics in the 
loop. There appears to be some correlation between the location of this 
peak and the low frequency complex conjugate poles of the transfer function 
relating torque exerted on the gimbal by the rotor to the input angular rate 
(Eq. (4. 1-7)). The peak gain is centered at about 175 Hz, and the poles of 
the transfer function are at about 190 Hz. From a practical viewpoint this 
change in the frequency response can constitute a serious limitation on the 
useful bandwidth of the gyro loop. In the gyro analyzed, greater rotor bear- 
ing rigidity is required if the peak is to be moved to a higher frequency and a 
larger bandwidth (predicted by assuming simple gyro dynamics) is to be 
realized. 
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Figure 5. 2-5 Frequency Response of a Ternary Gyro: 

Comparison of Analytic and Simulation 
Results when Rotor Bearing Compliance 
is Considered 
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5.3 EFFECT OF RANDOM SIGNAL AT THE NONLINEARITY INPUT 

Describing function analysis is used to obtain the closed loop 
response of a ternary gyro when a Gaussian random signal is present at 
the nonlinearity input along with a sinusoidal signal. The random signal, 
described by its standard deviation, o r , can be present as a result of a 
random input to the gyro, noise in the loop, or deliberate injection in the 
form of a random dither. The present discussion makes no assumption 
about the origin of the random signal. 

The Describing Function - Due to the presence of a random signal 
in addition to the sinusoid, the describing function gain to a sinusoid is dif- 
ferent from that used in the previous section. It depends not only on the 
amplitude A of the sinusoidal input to the nonlinearity but also on the magni- 
tude of the random input. A random signal can be described by its statistical 
moments. If it is assumed to be a Gaussian, zero mean variable its distri- 
bution is completely specified by a standard deviation, a r „ The describing 
function for the sinusoid is written as N^(A,a r ) and is given for the ternary 
nonlinearity in Appendix E of Ref. 3. As o r approaches zero N^(A, o r ) be- 
comes the describing function for a pure sinusoidal input. Figure 5. 3-1 
illustrates the effect of the random input on the sinusoid transfer charac- 
teristics of the ternary nonlinearity. It is similar to Fig. 5. 1-4. The out- 
put characteristic shown is the average amplitude of the output harmonic 
of the same frequency as the sinusoidal input. The effect of a random input 
to the ternary nonlinearity is to make the rise of the sinusoidal input-output 
curve more gradual. When a r >0 and A<6 there is a finite probability 
that the net input to the nonlinearity will exceed 6 and have the same sign 
as the sinusoid. There is a smaller probability that it will exceed the dead- 
zone with the opposite sign. When the mean is taken, a non- zero average 
value for the appropriate nonlinearity output harmonic results. In this way 
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the sinusoid is boosted through the deadzone by the random signal. On the 
other hand, when A exceeds 6 the fact that the noise input to the nonlinearity 
can have a sign opposite to that of the sinusoid produces a describing func- 
tion gain somewhat lower than that occurring in the absence of random sig- 
nals. It can be seen from Fig. 5. 3-1 that the net effect of adding a random 
signal at the nonlinearity input is to make the describing function gain for 
the sinusoid less amplitude- dependent (more nearly linear). 



Figure 5.3-1 Effect of a Random Input on the Sinusoid 

Transmission Characteristics of a 
Ternary Nonlinearity 

Unfortunately, the presence of the random signal in the loop 
involves more than just a change in N^. Strictly speaking, c r r must be 
determined through the use of a describing function gain N r (A,u r ) and 
knowledge of the characteristics of the random signal at the point where 
it enters the gyro loop. Since N A and N r both depend on A and cr r , solu- 
tion of the closed loop relations in the presence of both a sinusoid and a 
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random signal requires a complicated iterative procedure. However, in 
order to illustrate the effect of random inputs to the nonlinearity we will 
assume different values of a r without specifying their origin. 

The fact that the curves obtained represent the average response 
to a given sinusoid does not detract from their value. If we are interested 
in the effect of gyro output on a resonant mode of the direction cosine ma- 
trix calculation or on the system drift rate generated, the mean value 
representation is often as useful as an exact one. 



Figure 5.3-2 Ellipse Intersection Points Caused 

by Random Nonlinearity Input 

Changes in the Frequency Response - Figure 5.3-2 indicates 
the ellipse intersection points for the sinusoid input-output curve which 
result when a random signal also enters the nonlinearity. Note that the 


97 


likelihood of multiple intersection points occurring on any one ellipse is 
reduced. In this particular example there are no multiple intersection 
points. Also, because there can be an output for A <6, all ellipses have 
at least one intersection point. The effect on the sinusoid frequency re- 
sponse is illustrated in Fig. 5.3-3. Since no ellipse has multiple inter- 
sections, there is only one value of closed loop gain at each frequency; 
the jump resonance phenomenon is eliminated. Also , the closed loop gain 
is more than zero for all frequencies and the sharp cutoff found in Fig. 

5. 1-4 is eliminated. Figure 5.3-4 shows the three frequency responses 
corresponding to the three nonlinearity curves of Fig. 5.3-1. As cr r in- 
creases, the jump resonance disappears and bandwidth increases. Gradual 
gain reduction at high frequencies replaces the sharp attenuation. In short, 
when the random signal is present, the closed loop frequency response 
appears to be more like that of a linear system. 


98 






6 . 


GYRO AND SYSTEM ERROR SIMULATION 


SUMMARY 


A hybrid simulation which will facilitate the study of 
errors in strapdown gyros and strapdown navigation 
systems is described. To begin with, the motivation 
for a strapdown gyro and error simulation is discussed. 
The form of the proposed simulation is presented, along 
with some reasons for rejecting other approaches to em- 
ulating strapdown gyro and system behavior. The section 
closes with a description of that portion of the hybrid 
simulation which computes a measure of system attitude 
error from errors at the outputs of the individual gyros. 


6. 1 MOTIVATION FOR A SIMULATION 

The analysis and results reported in Ref. 1 and 2 are based on 
several assumptions and simplifications. Before proceeding with further 
work based on those results, it would be prudent to establish some check 
on the validity of these approximations. While simulation is not a perfect 
substitute for the construction and testing of hardware it is considerably 
less expensive and far more versatile. As an intermediate step between 
analysis and prototype construction, simulation is ideally suited to the 
present stage of the investigation. The value of simulation has already 
been well demonstrated in this report- -on the basis of the comparisons 
obtained in Sections 4 and 5, the response characteristics of pulse rebal- 
anced gyros using delta modulation torquing can be determined analytically 
with a considerable amount of confidence. 
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Other topics which were treated in an approximate manner in 
the above references include crosscoupling errors and random vehicle 
angular motions. When crosscoupling errors were computed in the case 
of oscillatory angular motion, the gimbal output angle was obtained from 
a linearized model for the binary pulse torqued gyro. The linear sensor 
approximation does not take into account the pulsed nature of the rebalance 
torque or float angles generated by limit cycling. In addition, treatment 
of crosscoupling errors in the ternary pulse torqued gyro is very difficult 
even with the aid of describing function techniques. Because of the non- 
linear nature of pulse rebalanced gyros , no confidence can be placed in 
the analysis of crosscoupling effects unless the results are substantiated 
by simulation. 

In the portion of Ref. 1 dealing with gyro parameter optimiza- 
tion, several approximations and simplifications were made to permit 
easy calculation of gyro errors caused by random vehicle angular motions. 
Essentially, these assumptions involved the use of sinusoid transfer 
characteristics predicted by the analysis of Section 4 and 5 to compute 
errors induced by random vibrations. A simulation which generates real- 
istic analogs of random angular motion about all major sensor axes is 
needed to check the assumptions. 

Several questions which arose and were left generally unanswered 
in the preceding work can also be treated easily by a well designed simula- 
tion. In the discussion of strapdown gyro error compensation, the possi- 
bility of using the three basic gyros to measure angular motion for their 
own compensation is mentioned. Stability and accuracy analyses of this 
approach are very difficult to perform analytically because many closed 
loops, often containing time- varying gains, are generated. However, the 
technique is so attractive in terms of cost, power, weight and reliability 
that it must not be ignored. A well designed simulation of three gyros and 
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the compensation calculation will permit an empirical determination of the 
feasibility of this approach. In addition, simulation provides an inexpensive 
way of checking other proposed gyro error compensation schemes, includ- 
ing their sensitivity to the accuracy of various motion sensors. 

In the work reported in Refs. 1 and 2, a system level figure-of- 
merit was developed to indicate the effect of gyro errors on strapdown 
system accuracy. Equations were derived relating the time histories of 
gyro errors to system attitude error. These expressions were exercised 
in the case of certain easily described gyro output errors (bias, white noise, 
oscillations, etc. ) to give simple analytic relations between sensor and 
system drift rates. Whenever examples were presented and parameter 
optimization was illustrated, the latter relations were used rather than the 
original, more complex equations. The approximations which resulted 
should be substantiated by simulation. 

In addition to verifying results from prior analyses, a well 
designed flexible simulation will be invaluable in the continuing study of 
strapdown sensors. Compensation within the gyro loop is easily analyzed 
and the effects of various torque generator characteristics on system 
errors can be studied. Test procedures for determining basic sensor 
error parameters can be evaluated. Using the simulation, dynamic test 
inputs can be designed and simulated information will be available for use 
in comparing different test data reduction techniques. Finally, the versa- 
tility and flexibility of a simulation must be emphasized. While the ulti- 
mate verification still lies in construction and testing of actual hardware, 
when it is desirable to investigate many proposed changes in gyro para- 
meters only simulation provides an inexpensive and expedient means. 
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Description of the Projected Simulation - The simulation planned 
and under construction is illustrated in Fig. 6. 1-1. Investigation indicates 
that a purely digital approach will be too slow if numerical integration is 
used to simulate the complex, high frequency dynamics to be considered. 

On the other hand, system error is to be computed using the error growth 
difference equations described in Ref. 2. This requires a digital computer. 
The hybrid simulation shown permits simultaneous determination of high 
frequency variables and accurate calculation of system error growth. 

Three gyros, including pulse rebalance torquing electronics, will be simu- 
lated on the analog portion along with random and deterministic motion in 
three dimensions, gyro motion- induced error torques and error compen- 
sation schemes. The gyro models will include structural dynamics when 
desired. Sensor outputs and true angular motion will be provided to the 
digital portion for calculation of system attitude errors. The simulation 
planned is a flexible tool for analyzing additional aspects of strapdown sen- 
sor errors. Gyro models can be changed to accommodate studies of exotic 
sensors and accelerometers can also be simulated. 


ANALOG COMPUTER 


r 



J 


Figure 6. 1-1 Hybrid Simulation of Gyros and Calculation 

of System Attitude Error 
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6.2 PROGRAM FOR THE DIGITAL PORTION 


One function of an inertial navigation system is to provide an 
indication of vehicle orientation with respect to a reference coordinate 
frame. In strapdown systems this requirement is satisfied by computing a 
direction cosine matrix, C, relating vehicle and navigation reference 
frames. Given an initial C matrix, the direction cosines can be computed 
from the history of rotations about each body axis. If we compute C at 
discrete intervals T seconds apart, the direction cosines at the end of the 
(n + l)st interval can be expressed in terms of those at the end of the prior 
period and a matrix of incremental angles, A0 (see Ref. 1). The matrix 
A0 contains the angles through which the body axes have rotated since the 
last sample time. The equation relating the two direction cosine matrices 
is: 


'n + 1 




( 6 . 2 - 1 ) 


when computer considerations such as quantization, word length and trun 
cation are ignored. A0 n is defined by 


/*(n + 1)T 

0 



> 

II 

w z 

0 

-C0 X 

nT 

_" w y 

«x 

0 


( 6 . 2 - 2 ) 


where o)y and a> z are the angular rates about the body axes. In prac- 
tice the elements of A0n are measured by the x, y and z-axis gyros. 

Generally, there is a difference between the actual A0 n matrix 
and that obtained from gyro measurements. The difference between the 
actual A0 n matrix and that given by the gyros is a third matrix, E n . The 
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presence of E n results in an error, AC, in the C matrix computed from 
the gyro measurements. Given E n for all values of n, we can compute AC 
at the end of each interval. A difference equation for AC which retains 
second order terms in small matrices is obtained from the series expan- 
sion of Eq. (6.2-1): 

AC n+i = C n+1 (E n+ l(E n ) 2 ) + AC n C*C n+1 (I + E n ) (6.2-3) 

A figure- of- merit for gyro performance, p n , can be computed from AC n 
(see Ref. 1): 

p n = \ trace (AcJ AC n ) (6. 2-4) 


Description of the Program - At each sample time the digital 
computer portion of the simulation receives the simulated gyro outputs 
(pulse counts) and true body rotations (actual A0 n matrix) as inputs. The 
program then proceeds to compute E n , C n + j, AC n+ \ and p n . The ma- 
trices C n + i and AC n + j are stored for the next cycle and p n is provided 
as an output. Figure 6. 2-1 is a flowchart for the program. On this chart, 
C Q and AC 0 are initial condition matrices for C and AC, respectively. 

The number n is the current number of compute cycles executed. This 
program has been run on a digital computer identical to the one to be used 
in the simulation. The compute cycle was found to be about 45 msec in 
duration. This constitutes a lower limit on the size of the sample period 
T if the digital computer is to process data when received and the analog 
simulation is not time scaled. The program was also checked out on an 
IBM QUIKTRAN time- sharing system. Simulated gyro error inputs in the 
form of ramps, sinusoids and coning motion were used. Jn every case the 
computer program output (p n ) agreed with the expected result. The program 
is considered ready for use. 
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STOP 

Figure 6. 2-1 Flow Chart for the Digital Part 

of the Hybrid Simulation 
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7. 


CONCLUSION 


7.1 CONCLUSIONS 

Progress has been made in developing simple and accurate 
analytic descriptions of strapdown gyro behavior. A gyro employing bi- 
nary delta modulation torquing was represented by a linear transfer func- 
tion and the approximation was verified by simulation. The only time that 
a binary gyro is not well represented by a linear transfer function is when 
the limit cycle is extinguished by an input of sufficient frequency and amp- 
litude. However, this situation is easily predicted by analysis and inputs 
required to quench the limit cycle are in any case unlikely. It was shown 
that gyros incorporating the time modulation rebalance scheme can also 
be represented by a linear transfer function. 

Excellent agreement was found between the transfer character- 
istics predicted analytically for a ternary pulse rebalanced gyro and those 
experienced by simulation. The analysis in the case of the ternary gyro 
is more complex than for the other two sensors, reflecting in part the 
nonlinear nature of this gyro loop. The amplitude dependent nature of the 
ternary gyro response is well illustrated in Section 5. 

The effects of complex dynamics within the gyro loop were ac- 
curately predicted analytically. The describing function techniques used 
throughout are capable of handling linear loop elements of any complexity. 
Consequently, including a linear model of rotor bearing compliance in the 
formulation of the analysis did not present any great difficulty. It was 
shown both analytically and by simulation that rotor bearing compliance 
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can alter gyro transfer characteristics significantly. However, it is still 
possible in the case of the binary gyro to represent the sensor response 
to sinusoidal inputs by a linear transfer function. Sampling within the gyro 
loop was also shown to reduce the instrument’s bandwidth. 

Most of the approximations regarding gyro loop behavior discussed 
in this report pertain to sinusoidal inputs only. To analyze the contribution 
of random vibrational motion to strapdown system errors it is necessary 
to describe the transmission characteristics of gyro loops when random 
inputs are applied. For gyros using binary delta modulated and time modu- 
lated torquing the answer seems to be a simple extension of the work 
reported above. Since the sinusoidal response of these two sensors is 
linear the same loop descriptions may suffice. The problem of developing 
a simple expression describing the propagation of random signals through 
a ternary gyro appears to be much more difficult; the effort required to 
investigate response to sinusoidal inputs is considerable and there is no 
reason to expect that random inputs will simplify the problem. The ques- 
tion of higher harmonics in the gyro output also remains; recall that the 
definition of frequency response involved only that component of gyro output 
which had the same frequency as the input. The problem of higher har- 
monics present in the square wave form of the gyro output must eventually 
be treated. In any case, the presence of these extra outputs does not in 
any way detract from the importance and validity of the excellent prediction 
of loop transfer characteristics demonstrated. 

Finally, all of the results obtained this far relate to the gyro 
output. In order to compute crosscoupling errors the float angle response 
must be described. This will require additional, though probably not 
radically different, approximations to the behavior of nonlinear gyro loop. 
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7. 2 CONTINUATION OF EFFORT 


Development of the gyro error simulation will be the focus of 
future work. The time modulation gyro torquing scheme will be simulated 
and frequency response characteristics obtained will be compared with 
those predicted analytically. Consequently, it should be possible to verify 
the analytic description of this kind of gyro loop presented in Section 4 of 
this report. Analogs of motion-induced error torques will be applied to 
the gyro models to study the propagation of errors in the individual gyros 
and in a strapdown system. The sinusoidal-input transfer characteristics 
discussed in this report will be used to predict transmission of random 
signals through the gyros and simulation will be used to check the analysis. 
In this manner, it is expected that the linear approximations to the binary 
and time modulation torqued gyros will be certified for use in analysis of 
random vibration- induced errors. 

The possibility of compensating environment- induced gyro 
errors by processing outputs from the basic sensor triad will be explored. 
Emphasis will be placed on compensation through corrective torquing but 
some effort will be devoted to the problem of compensating at the gyro 
output. It is anticipated that the accuracy potential of both these techniques 
will be determined by simulation in the near future. If certain aspects of 
gyro error compensation (for example, removing output axis angular 
acceleration errors) are shown to be unsatisfactory when only the gyro 
outputs are processed, the simulation study will be extended to cover the 
use of additional sensors. 

Future analytic work will be directed toward describing the 
propagation of dynamic errors through a set of ternary pulse torqued 
gyros and into system misalignment. In particular, it is thought that 
crosscoupling errors in the ternary gyro may be treated by the 
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describing function approach. In brief, the future effort is expected to 
advance understanding of strapdown sensor errors and their reduction 
toward the point where it can be used as the basis for some hardware 
construction decisions. 



APPENDIX A 


GIMBAL MISALIGNMENTS 
CAUSED BY OUTPUT AXIS ANGULAR MOTION 


A single- degree- of- freedom gyro experiences torques about its 
input axis which result from angular rates about the output axis. These 
torques misalign the float about the input axis, IA, to a degree which de- 
pends on the stiffness of the float suspension. Suspension stiffness to- 
gether with the floatation fluid determine the time constant of the relative 
motion between case and float. In this section we analyze the time response 
of the float misalignment about IA to a step angular rate applied about the 
output axis, OA. 


SRA 



o 


Figure A-l Model for Gyro Case and Float Dynamics 
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Figure A- 1 shows the model assumed for float dynamics about 
the input axis. The suspension stiffness is represented by linear springs 
and fluid damping by a pair of dashpots. 

q = Float- to- case angular misalignment 
about IA 

ot 0 = Float- to- case angular misalignment 
about OA 

coj = Case angular rate about IA 

co 0 = Case angular rate about OA 

lii = Float moment of inertia about IA 

H = Wheel angular momentum 

A = Distance between float center of mass 
and each float bearing in cm 

D = Equivalent linear damping constant in 
dyne-sec/cm 

Cj = Equivalent rotational damping constant 
in dyne- cm- sec 

= DA 2 

K = Equivalent linear spring constant in dyne/cm 

The equation for angular motion of the float about IA is derived: 

Torque on Float = Gyroscopic Torque - (Spring &Damper Reaction Torques) 
= L. (wj + a.) 
or 

L. (w. + o^) = -H (w 0 + a Q ) - 2Da\ - 2KA 2 a . (A-l) 

Simplifying, we have 

- + a Q ) - w. = + 2 CGjjOj + 0 2 Oj (A-2) 
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where 



and 



If we assume co. and a Q to be zero, Eq. (A-2) reduces to 


V*cVi 

11 


(A-3) 


(A-4) 


(A- 5) 


For most gyros the float motion is heavily damped and the 
response to a step of angular rate about OA can be accurately described by 
a single time constant, t' 



(A-6) 


r = 


1 

n „(t -VC 2 ^) 



(A- 7) 
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Equations (A-3), (A-4), (A- 6) and (A- 7) are used to compute float misalign- 
ment behavior for different values of Cj and K. The fixed gyro parameters 
chosen are 


H = 2 0 5 x 10^ dyne-cm-sec 

4 - 3.6 cm 

Ijj = 300 gm-cm^ 

Three plots are presented to assist the gyro designer in picking the appro- 
priate suspension stiffness and damping coefficient to constrain misalign- 
ment about IAo 

Figure A- 2 shows the behavior of the float time constant, r, as 

a function of K and Cj. Given r we can enter Fig. A-3 and obtain a time 

€% 

history of the normalized float angle, 24 Ka./Ha> 0o Note that any time 
constant obtained from Fig. A- 2 implies 4 = 3. 6 cm. Figure A-4 allows 
us to obtain the misalignment angle as a function of time directly. The 
effect of suspension stiffness on the time response is readily observed. 

If the suspension is not stiff enough for the range of expected 
output axis rates, the float pivots can strike the bearings. The resulting 
static friction between pivot and bearing surfaces provides a static break- 
out torque, t s , which masks small torques about the output axis caused by 
input angular rates. Accordingly, when the three- level torquer is used 
input rates smaller than r s /H are not seen at the sensor loop output, and 
the gyro has an effective rate deadzone. The force, F, pushing each pivot 
against its bearing results from the gyroscopic torque caused by an output 
axis angular rate and from forces applied by the gimbal suspension: 



K (dyne/cm) 


Figure A- 2 


Float Misalignment Time Constant as a 
Function of Suspension Stiffness K and 
Damping Constant Cj 










K (dyne/cm) 

Figure A-4 Normalized Misalignment for Values 
of Suspension Stiffness 
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where d is the clearance between bearing and pivot. Knowing M static > 9X1 
effective measure of static friction at the pivot-bearing interface, we can 
evaluate r s : 


r s = ^static F 

/ Hw o \ 

- ^static " 2Kdj 

The rate dead-zone corresponding to this break-out torque is 

/ w o 2Kd\ 

± ^static \ i ~ h / 

In a binary gyro the constant application of torque by the torque generator 
prevents this dead- zone. 
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APPENDIX B 


MOTION INDUCED ERRORS 
IN THE SINGLE-DEGREE-OF-FREEDOM 
PENDULOUS ACCELEROMETER 


Errors introduced in the single -degree-of-freedom (SDF) pen- 
dulous accelerometer are discussed below. Whenever possible the error 
torque mechanisms are related to those given for the single-degree-of- 
freedom gyroscope in Ref. 1. Free use is made of both the derivations and 
discussions presented in the reference. 

The SDF pendulous accelerometer is illustrated in Fig. B-l. 

Two major differences between this representation of the instrument and 
that presented for the SDF gyro are obvious. The direction perpendicular 
to the output and input directions is called the pendulum (p) axis rather than 
the spin (s) axis. Also, the instrument is assumed to consist of only two 
basic parts: a case and a combination gimbal and pendulum. The error 
torques induced in the SDF pendulous accelerometer by angular motion can 
be expressed by modifying the equations developed in Appendix A of Ref. 1. 
As a result of the instrument model shown in Fig. B-l, all terms in 
Eq. (A- 19) of the reference which involve # Q and 0. are dropped, and the 
subscript g is deleted because the gimbal and pendulum are considered as a 
single rigid body. The subscript p is substituted for s in the remaining 
terms. Also, since there is no spinning wheel involved in the accelerometer 
under consideration, all terms involving angular momentum, H, and the spin 
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Figure B-l Exploded View of Single- Degree-of- Freedom 
Pendulous Accelerometer Showing Small 
Misalignment Angle, a 
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rate of the rotor relative to the gimbal, £2 , are deleted. The resulting 

s 

expression for output axis torques produced by angular motion is; 


W = I (a + 6d^+(l.. - I ^60 60. + I (oi - 60. co. ^ 

, 00 V O 0/ \ 11 pp/ pi op\ P 0 1/ 


'oifii* u o“p) + 1 pi{% - *i) 
oc ( I - I..)(co 2 - co. 2N ) 

0\pp 11/ \ p 1 / 


- CL (l - I.. ICO CO + I 60. 
^)|\pp 11/ op 00 lj 

+ CL. ( I - I.. ICO CO. + I CO 

Vpp 11/ 0 1 oo pj 


(B-l) 


Since the ideal accelerometer is insensitive to angular motion, all 
of the terms in Eq. (B-l) must be considered as error torques in the SDF 
pendulous accelerometer. However, as in the case of the SDF gyro, the first 
term, I qq a Q , together with any damping torque about the instrument output 
axis, is usually considered as part of the unavoidable sensor dynamics and 
included in any ’’ideal" SDF pendulous accelerometer model. The remaining 
error terms can be divided into several broad categories similar to many 
exhibited by the gyro. Sensitivity to angular accelerations is present. The 
principal contribution, that caused by angular acceleration about the sensor 
output axis, is unavoidable because of the nature of the pendulous acceleration 
sensing instrument. Several anisoinertia terms also appear. Because the 
gimbal and pendulum of the accelerometer can be viewed as a rigid body, the 
opportunity exists to design an instrument which experiences no error torques 
from this cause. Finally, product of inertia terms appear in Eq. (B-l). The 
first of these terms can be larger than that experienced by a SDF gyro if the 
pendulum shown in Fig. B-l deviates from the p' direction in the o'-p' plane. 


123 



Linear Motion - Error torques generated by linear motion 
(acceleration) of the SDF pendulous accelerometer can be established using 
the error torque equation* 

[M f ] , = m(6 x [a] Q ) • o’ (B-2) 

where m is the mass of the gimbal and pendulum combination, 6 is the vector 
of displacement between the gimbal-pendulum center of mass and the output 
axis expressed in gimbal coordinates, and [a]^, is the linear acceleration 
vector, similarly expressed. Notation will be as near to that of Section 2 
and Appendix A of Ref.l as possible. Acceleration expressed in the gimbal 
frame is related to that experienced in the case frame by: 


[a ) Q = [a] c - a x [a) c (B-3) 

Substituting Eq. (B-3) into Eq. (B-2), the output axis torque caused by the 
linear acceleration of a pendulous accelerometer is: 



a. + a a 
1 op 


- Vo) " m6 i( a p + “i a o ■ Vi) 


(B-4) 


The effect of linear compliance on accelerometer errors can be 
illustrated by describing the center of mass displacement as a function of 
linear acceleration: 


6 = m[K .a. +K a +K a |+ 
p V pi i po o pp VI 

6. = m(K..a. +K. a +K. a )+ & 
i \ 11 1 10 ° ip P f o 
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Substituting Eq. (B-5) into Eq. (B-4), a more detailed torque equation 
results: 

[M f ] § = m6 po (a* + 0 ^- 0 ^) 

- (ap + cx^Q- a^) (B-6) 

+ m ^pi 8 ! + ^po^i a o + ^pp ~ ^ii^ a i‘*p ” ^io^o^p ~ ^ip^J 


The first term of Eq. (B-6), m6p a| , measures linear accelera- 
tion along the input axis. This is the only output axis torque in the ideal 

pendulous accelerometer. The pendulosity m6„ is designed into the instru- 

*o 

ment with care. All the remaining terms in this equation contribute errors 
to the accelerometer. The term m6p a 0 ap is basically a cross- coupling 
error arising from rotation about the single axis of freedom and m6 p a p a 0 
results from gimbal-to-case misalignment. Since accelerations along the 
input axis will cause considerable excursions of the gimbal angle, a Q , from 
null, sizeable rectification errors can be produced in this instrument by 
properly phased linear vibrations with components along the input and 
pendulum axes. The second line of Eq. (B-6) illustrates error torque 
contributions from unwanted mass unbalance and the last line expresses 
compliance error terms. It can be seen that linear compliance effects 
can produce constant error torques. The accelerometer is also subject 
to random error torques similar to those in an SDF gyro. The error in 
indicating linear accelerations along the case fixed input axis of an SDF 
pendulous accelerometer is simply the sum of all error torques, divided 
by the pendulosity, m6p Q . 
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Rebalance Loop Errors - Errors in the single- degree-of- 
freedom pendulous accelerometer which result from the rebalance torquing 
mechanism are identical with those described for the SDF gyro in Section 
2.4 of Ref. 1. The rebalance loop configurations commonly used are 
similar to those for the gyro and the input to the rebalance loop of either 
instrument can be viewed as a torque instead of an angular rate or a linear 
acceleration. 
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APPENDIX C 


ACCURACY OF THE SMALL SIGNAL APPROXIMATION TO THE 
TWO-SINUSOID-INPUT DESCRIBING FUNCTION 
FOR THE BINARY NONLINEARITY 


In Section 4. 1 the Two-Sinusoid-Input Describing Function (TSIDF) 
gain of the binary nonlinearity is used. The describing function gain for a 
small signal sinusoid in the presence of a large limit cycle sinusoid is ex- 
pressed there by the approximation: 

jj sj ^ (C-1) 

N B nA 1} 

where A is the limit cycle amplitude at the nonlinearity input and D is the 
relay drive (output) level. The describing function gain for the limit cycle 
in the presence of the small signal is given by 

n a* IS • (c - 2) 

It is desirable to evaluate the accuracy of each approximation. 

The required TSIDF’s are given in Appendix D of Ref. 3. De- 
noting the amplitude of the small signal as B, where B/A = p f p < 1: 


n b (a, B) = -|— 

ir Bp 


E (p) - (l - p 2 )K(p) 


(C-3) 


127 



and 


N a (A, B) = -^-E (p) 
it A 


(C-4) 


where 


K(p) 



dip 

2 . 2 , 
- p sin 0 


(Elliptic Integral of the First Kind) 


E(p) 



2 2 

- p sin ip dip 


(Elliptic Integral of the 
Second Kind) 


The percentage errors in TS3DF estimates caused by using the appropriate 
approximations are: 


% Error in N fi = 


Approximate Value of N R 


Exact Value of N 


- 1 


B 


x 100 


= 25 


i rp 2 - 4 [e (p) - (l - p 2 )K(p)] 
E (p) - (l - p 2 )K(p) 


(C-5) 


and 


Error in N A = 
A 


Approximate Value of 
Exact Value of N . 


- 1 


x 100 


= 50 


it - 2 E (p) 

— eTp) ' 


(C-6) 
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Percent Error in N^, 


Note that these errors are function of p only. Figure C-l shows the 
errors plotted as a function of p. As expected, the errors in the approxi- 
mations vanish as the amplitude ratio p goes to zero. 


p = ^ where: B = Signal Amplitude Entering 

A the Nonlinearity 

A = Limit Cycle Amplitude 

Entering the Nonlinearity 



Figure C- 1 Error in the Approximation 
to the TSDDF 
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APPENDIX D 


LIMIT CYCLE QUENCHING 
IN THE BINARY GYRO LOOP 


When determining the limit cycle of a binary pulse rebalanced 
gyro, the nonlinearity can be replaced by an equivalent describing function 
gain, N^o If both a signal, B sin vt, and the limit cycle oscillation, 

A sin out, are present at the input to the nonlinearity, Na remains es- 
sentially independent of the signal as long as B << A (see Appendix C). 

In this appendix the influence of a large signal on the describing function 
gain to the limit cycle is considered. The material presented follows 
Ref. 3, Example 5.3-2. 

If a limit cycle is to exist in a binary gyro, the following amp- 
plitude condition must be satisfied: 

Na I H(j«i) I = 1 (D-l) 

where H(s) is the open loop transfer function of the linear part of the gyro 
loop. The limit cycle frequency cty is determined from the phase condition: 

/H(jcoe) = -180° (D-2) 

Several possible values for <«$, may exist, but they are fixed for any given 
gyro. Consequently a limit cycle cannot exist if 

Na S I H(jw/) | t 15 - 3 ) 

for all possible values of cog . 
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The describing function gain, N^, of the binary nonlinearity is 
in general dependent on A and B. From the theory of Two- Sinusoid- Input 
Describing Functions (TSIDF’s), 

Na(A,B) = E(|) (D-4) 

IT 2 A A 

where D is the nonlinearity drive level, A/B > 1 and E(B/A) is the complete 

elliptic integral of the first kind (see Ref. 3). Maximizing over all 

values of A yields Na as a function of B: 

J A max 


n a . 


0.855D 


max 


B 

A 


= 0.909 


(D-5) 


Combining Eqs. (D-l) and (D-5), the value of the signal amplitude which 
causes the limit cycle to be ’’quenched,” B q , is found: 

B q = 0.855D | H(jcc £ ) | (D-6) 

The nonlinearity describing function gain to the signal when the limit cycle 
is quenched is found by substituting B/A = 0. 909 into the relation for 
Nj$(A,B) given in Appendix C and employing Eq. (D-6): 


N b (A,B)i _ 2.089D 

! B = B q if B q 

0.778 
I H(jWi) I 


(D-7) 


The binary torqued gyro loop can now be viewed as the system shown in 
Fig. D-l. 
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or 0 = B sin vt < A sin 



Figure D- 1 Binary Loop Viewed at the Point where 

the Signal Quenches the Limit Cycle 


In order to determine the values of input frequency and ampli- 
tude which cause the limit cycle to be quenched, the transfer function 
between input angular rate and float angle ot Q is written in terms of the 
transfer function for the various blocks in Fig. D-l. 
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From Eqs. (D-6), (D-7) and (D-8), the solution for the input amplitude 
required to quench the limit cycle can be written: 


^quench ~ 




a, 


(s) Bq 
s = \v 


1 + 


M. 


“°(s)>(s) 0.778 


net ‘te lH(i^)l 


M l, , “o 

10 . ^ M . 
i net 


(s) 


(D-9) 


0.855D|H(ioJjj) | 


s=]v 


The most convenient way to display Eq. (D-9) is a plot of ^q uenc h as a 
function of v. For a particular set of parameters, a single- valued curve 
results, delineating the boundary between the regions where quenching 
does and does not occur. Figure D-2 indicates how limit cycle quenching 
helps determine the theoretical frequency response curve seen in Fig. 

4. 2-4 for the binary loop with simple gyro dynamics. The top half indicates 
the quenching curve for this system, showing the frequency where quench- 
ing occurs when a 0. 1 rad/sec oscillatory input rate is applied. The re- 
sulting frequency response with the effects of quenching shown in the lower half 
of the figure. When the limit cycle is quenched, the loop output is a square 
wave at the signal frequency with an amplitude equal to the nonlinearity 
drive level. Torquer current is equivalent to 1 rad/sec. The fundamental 
component of the output has an amplitude of 4 /it = 1.276 rads/sec and the 
gain of the loop is (1. 276/0. 1) = 12. 76. The effect of quenching is illus- 
trated in Fig. D-3 for the loop with compliance dynamics discussed in 
Section 4. 3. 2. 
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Frequency (Hz) 


Figure D-2 


Limit Cycle Quenching and Its Effect 
on Frequency Response (Second 
Order Gyro Dynamics) 






Frequency (Hz I 


Figure D-3 Limit Cycle Quenching and Its Effect 
on Frequency Response (Rotor Bear- 
ing Compliance Considered) 
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Figure D-4 presents simulation records illustrating limit cycle 
quenching. The situation shown corresponds to suppression of the limit 
cycle as input frequency is increased to about 140 Hz, indicated in Fig. 
D-3. The frequency of wj was increased slowly but the change is not 
perceptible in the figure. Limit cycle quenching is seen as an abrupt 
change in loop behavior as reflected in the float angle. 


100 msec 


0. 1 rad/sec 

L_ 


n 


(a) Input Angular Rate, 


~r 

1 mrad 

_L 





(b) Float Angle, a Q 

Figure D-4 Binary Gyro Simulation Record 
Showing Limit Cycle Quenching 
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APPENDIX E 


DERIVATION OF THE ELLIPSES USED IN DETERMINING 
THE FREQUENCY RESPONSE OF THE 
TERNARY GYRO LOOP 


In analysis of the ternary pulse rebalanced gyro loop (see Sect. 5), 
a graphical solution is used. This is done by using two separate descriptions 
for the amplitude ratio between sinusoids on either side of the nonlinear element. 
One description, developed from linear analysis techniques, generates a 
family of ellipses. The equation for these ellipses is derived below. 


Figure E-l shows a generalized loop containing a single nonlinearity. 
N is an equivalent sinusoidal gain (describing function gain) for the nonlinearity 
and is assumed to involve no phase shift. The loop is driven by a sinusoid of 
constant amplitude, M, and frequency, co. The amplitude of the transfer 
function between the signals A sin (cot + <j>) and M sin cot is: 


Using the identity 


A_ 

M 


Pi 


J0 


1 + Np^/W 


e^ 0 = cose + jsinG 


Eq. (E -1) may be written as: 


(E-l) 


(E-2) 



Pi 


^ (l+Np 1 p 2 cos(0 1 +0 2 ))^ + N 2 Pi p 2 sin^(0 1 +0 2 )' 


(E-3) 
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M sin u>t -rQ 


AN sin(u>t t <b) 



L 


1 


L 


2 




N = equivalent sinusoidal gain 
of the nonlinearity 


Figure E-l A Loop Containing One Nonlinear Element 


Squaring and cross- multiplying gives: 

A^l+NpjPgCOsfej+eg)) 2 + A 2 N^PjP2Sin 2 (ej+©2) = M 2 p 2 , (E-4) 


Solving for the loop output amplitude, AN, we obtain 


AN = 


-Acos(6 


p l p 2 




(E-5) 


In this expression M is given and 6 ^, Pg are determined from the 

frequency to. Thus, for fixed input amplitude and frequency we have a relation 
giving the amplitude of the nonlinearity output as a function of A, the amplitude 
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of its input. Note than AN is independent of the value of N. Plotting AN 
for values of A yields the ellipse shown in Fig. E-2. Critical points in 
the construction are illustrated. 


(-M/P2) 

(M/p 2 ) 


AN 



Figure E-2 Nonlinearity Input-Output Relation 

Derived from Analysis of the Linear 
Portion of the Loop Shown in Fig. E-l 
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APPENDIX F 


SIMULATION OF GYROS 

WITH A SIMPLE REPRESENTATION OF DYNAMICS 


This appendix describes simulation of binary and ternary pulse- 
torqued gyro loops where the single-degree-of-freedom gyros involved are 
characterized by second-order linear dynamic equations. All simulations 
were scaled to run at 1/100 real time. The basic signal flow diagram 
describing the gyros is given in Fig.F. 1-1. The simulations differ only in 
the form of the nonlinearity assumed. Appendix G describes the simulation 
of gyros with more complex dynamics. 


F. 1 BINARY GYRO SIMULATION 

When simulating the binary torqued loop with simple dynamics, 
certain blocks shown in Fig. 1. 1- 1 are represented by the following linear 
transfer functions: 


M^/ o>i = H dyne -cm- sec/rad 


“o /M net = ic ; s ' r f V - n “v/dyne-cm 


k te 

M a Ji L _ = — — , — t dyne - cm/ ma 


tg^tg sr tg + 1 
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Definitions: 

«i 

M net 

“o 

*tg 

Mtg 

Mi 


Input angular rate, in rad/sec 

Net torque about the gimbal output axis, in dyne- cm 

Signal generator output, in mV 

Torquer current input, in mA 

Torque r output, in dyne- cm 

Gyroscopic torque, in dyne- cm 


Figure F. 1-1 Signal Flow Representation of 

Gyro Loop Simulation 


144 








The analog mechanization of these relations is shown in Fig. 

F. 1-2. Pot settings are given as algebraic combinations of loop parameters. 
To facilitate interpretation of analog voltages as physical quantities, a list 
of scaling factors is presented. The binary nonlinearity mechanization is 
given in Fig. F. 1-3. 



cc. 0.1 rad/sec 

M net 1.66 x 10 4 dyne-cm 

a 0 1.66 mV 

i^g 1. 66 mA 

Note: D is the nonlinearity drive level in mA. 2 xlO 4 dyne-cm 

One second simulation time = 0.01 sec real time. 

Figure F. 1-2 Mechanization of the Linear Portion 

of the Gyro Loop 

Prior to measuring the binary loop frequency response, the 
system is allowed to limit cycle with no input, and the oscillation ampli- 
tude and frequency are compared with values predicted from describing 
function analysis. If limit cycle quenching (see Appendix D) occurs during 
frequency response measurements, the float angle waveform becomes 
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nearly sinusoidal and has a large amplitude. The output amplitude, as 
indicated by the Fourier Detector, rises to 4D /ir rads/sec, a value deter- 
mined only by loop parameters, and is insensitive to further frequency 
changes. 


Clock 



(a) Block Diagram Model (b) Simulation Mechanization 


Figure F. 1-3 Binary Nonlinearity Mechanization 

For any given amplitude, critical points in the region of limit 
cycle quenching are given by four numbers. These are the frequencies at 
which the limit cycle is observed to be quenched and to reappear as the 
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input frequency is raised, and the corresponding points which occur when 
the frequency is lowered. The four points determine a region of jump 
resonance on the frequency response curves. When the input frequency 
is well above that of the limit cycle, the output of the Fourier Detector 
is quite noisy and difficult to read. 


F. 2 TERNARY GYRO SIMULATION 

The basic difference between the binary and ternary rebalance 
loops is in the nonlinearity used; the linear transfer functions and the 
placement of elements within the loop are the same regardless of which 
torquing scheme is used. Accordingly, Fig. F. 1-2 also describes the 
linear part of the ternary loop. Mechanization of the ternary nonlinearity 
is shown in Fig. F.2-1. 

Before measuring the ternary loop frequency response, the 
discriminator level 6 is lowered until limit cycling starts. The maximum 
value of 6 causing oscillation and the oscillation frequency itself are com- 
pared with those obtained from describing function analysis. During the 
frequency response measurement, jump resonance appears as a sudden 
increase in float angle and Fourier Detector output for a small frequency 
change. As with the quenched binary loop, the gain changes associated 
with ternary jump resonance occur at frequencies which depend on the di- 
rection of approach. However, in the ternary gyro there is no region in 
which the gain is insensitive to frequency changes such as can occur during 
limit cycle quenching in the binary gyro. There is a frequency above which 
loop output ceases because the peak value of the float angle no longer ex- 
ceeds 6. Once this frequency is exceeded, the loop output remains at zero 
unless the input frequency is lowered to a certain value where the output 
becomes non- zero again. The frequency of this second point is known 
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analytically but can be difficult to check by simulation because computer 
noise causes the float angle to drift. To remove this source of error, the 
input frequency was increased until the gyro output disappeared. The ob- 
served float drift rate was then corrected. When this procedure is followed, 
the loop output behavior is very close to that predicted by analysis. 


Clock 



Switches 


Resistor 


(a) Block Diagram Model 


Networks 

(b) Simulation Mechanization 


Figure F.2-1 Ternary Nonlinearity Mechanization 
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APPENDIX G 


SIMULATION OF GYROS 
WITH ROTOR- TO-GIMBAL COMPLIANCE 


The gyro simulations which take into account rotor- to- gimbal 
compliance are illustrated by the flow diagram of Fig. F. 1-1. The tor- 
quer, nonlinearity and Fourier Analyzer portions of the simulations are 
identical with those described for simulating simpler models of pulse 
torqued gyros (see Figs. F. 1-2, F. 1-3 and F.2-1). Only the dynamic 
relations between Mj and to* and between a 0 and M ne j- are treated in 
detail below. Practical considerations in obtaining frequency response 
curves are also discussed. 


G. 1 SIMULATION MODEL INCLUDING ROTOR- TO- GIMBAL COMPLIANCE 


The transfer functions relating input angular rate, gimbal torque 
and gimbal angle, derived in Appendix I, are: 


Mj ^ - H (4p 2 y 2 s 2 + 4pi$s + 4) 

( s ) = ■ 3 2 , 2 2 2 V 2 . 3 4 

^ s +4pi> n s +(2^ n + 4p i^ n +03g)s +Apu n + v n 


(G. 1-1) 


a. 


^net 


(s) 


1 (s 4 + 4py n s 3 + \2v\ (2p 2 + 1) + u |]s 2 + 4py 3 s + V 4 ) 

sI ° 0 g (s 5 + C 4 s 4 + C3S 3 + C 2 s 2 + C 1 s + C 0 ) 

(G. 1-2) 
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where 


C 4- 4 ^ + l44 

a 2 2 2 4<w n k 2pU n d 

C 3 S 2^(2, 2 + l ) + M ^— +I _ + 1 _ 

c 4 4p( 3 “g 2^(2P 2 + D 2^ n k d 2 ^ 


n + r . 


r f 


"" u / A A\ 

^ + W (t ' n+ g) 


g 


°i - v * 


4 


0 


The analog computer mechanization diagram for Eq. (G. 1-1) appears in 
Fig. G. 1-1. Appropriate potentiometer settings are given as algebraic 
quantities. Figure G. 1-2 shows the mechanization diagram for Eq. (G. 1-2). 
Figure G. 1-3 is included to illustrate analog computer scaling of important 
signals in the simulation. 


G. 2 BINARY GYRO SIMULATION 

The binary loop with simple dynamics characteristically limits 
cycles at a single frequency because there is only one point at which the 
open loop phase shift is -180°. When rotor-to-gimbal compliance is con- 
sidered, there can be more modes, usually three, each one characterized 
by a different frequency and amplitude. It is shown in Section 4. 1 that the 
limit cycle frequency determines the gain of the nonlinearity to the signal. 
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Potentiometer Values: 




A 

H 


E = 

40, 3 
10 10 


166 x 10 4 


B 

2, 2 (2o 2 - 1) * 

2 

F = 

4o 2 ,; 

10 8 



10 5 


u* 



4 ^n 

c ■- 

n 


G = 



To 14 



10 10 

D = 

40 •'n 
10 3 


J - 

4 

v n 

10 14 


Figure G. 1-1 Mechanization Diagram for Eq. (G. 1-1) 
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G 

H 

I 

J 

K 

L 


K St! ' 

10' B 

<% loo,,' ■ 

10" 12 

( (21/ 2 (2 p 2 

1' • ' io" 6 

10 loo,. 


<< p ‘'n'Ioo ( ,l 

• 10' 7 

(4 QV n / Ioog i 

- 10" 1 


Figure G. 1-2 Mechanization of Eq„ (G. 1-2) 
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Consequently, if the limit cycle switches to a new mode, the loop charac- 
teristics as seen by the input signal can be radically altered. The loop 
1ms a different closed loop frequency response for each limit cycle mode 
encountered. 


1 volt = | 


1.66 mV (Binary) 
0.166 mV (Ternary) 



Figure G. 1-3 Computer Scaling Used in Simulation of 

Gyros with Rotor- to- Gimbal Compliance 


Before a frequency response simulation is made, all possible 
limit cycle modes are determined analytically. One check of the simula- 
tion is to verify the existence of the calculated limit cycle modes. In this 
manner information is obtained on how the loop must be stimulated in 
order to make the various modes appear. After the simulation has been 
verified as a true representation of the loop dynamics, a particular limit 
cycle mode is obtained and the frequency response for this mode is mea- 
sured. During the measurement care must be exercised to avoid those 
frequency regions where different modes may be excited. This procedure 
is repeated until all modes have been considered. 
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An additional characteristic appears, caused by the lightly damped 
complex poles of Eq. (G. 1-1), when rotor-to-gimbal compliance is considere 
As the input frequency nears the lower of the two pole frequencies, the input 
torque applied to the gyro loop increases until it reaches a level sufficient to 
quench the limit cycle. A peak in the frequency response occurs when this 
condition is reached. For the system dynamics considered to date, this 
response peak occurs at a lower frequency than any resonance attributable 
to the dynamics of the gyro torque loop. Since the usable frequency range 
of most gyro loops consists of the region below the first gain peak, there 
is little need, given the present dynamics , to examine the region above the 
first pole frequency. Thus, in the present case the important part of the 
closed loop frequency response is determined by rotor-to-gimbal compli- 
ance characteristics and not by the float and torquer time constants. 


G. 3 TERNARY SIMULATION CONSIDERING COMPLIANCE 

Most of the discussion in Appendix F concerning the ternary 
simulation applies when rotor-to-gimbal dynamics are considered. The 
discriminator level 6 is lowered in the simulation until the loop starts to 
limit cycle. The limit cycle mode(s) thus obtained are checked against 
those predicted by analysis. If the modes are as predicted the simulation 
represents the desired dynamics. As with the binary loop (Section G. 2), 
a pole of Eq. (G. 1-1) causes a peak which dominates the closed loop fre- 
quency response of the ternary loop. The only difference is that the peak 
is characterized by jump resonance rather than by quenching. 
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APPENDIX H 


THE FOURIER ANALYZER 


One problem encountered in measuring the sinusoid frequency 
response of the pulse rebalanced gyro loop is that the loop output is not a 
sinusoid. The loop output is taken at the output of the nonlinearity; it 
consists of a series of rectangular pulses. When an input of the form 
A sin tot is applied to the loop, the loop output y is of the form 


s c 

y = B 1 sin cot + B 1 cos cot + 


E 

n = 2 


(B s sinncot + B c cosncot 
\ n 


n 


). (H-l) 


There is no constant term at the output because the input is assumed to be 
an unbiased sinusoid and no rectification takes place within the loop. The 
loop transfer function evaluated at the frequency co is defined by the 
equation: 


Loop Transfer Function = 


Phasor Rep, of the Output Component at Freq. co 
Phasor Representation of the Input (Freq. co) 




(H-2) 
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It is necessary to measure the amplitudes B^ S and B^ c of the in-phase and 
the quadrature components at frequency to in order to calculate the closed 
loop transfer function at that point. 


H. 1 ANALYZER OPERATING PRINCIPLES 

The Fourier analyzer operating concept is based in part on the 
following two identities: 


sin 2 0 = ^ (1 

- cos 20) 

(H-3) 

cos 2 0= ^ (1 

+ cos 20) 

(H-4) 


If the gyro output is multiplied by A sin cot, the result is as follows: 


( S C ' 

B n A sin tot sin ncot + B n A sin tot cos ntot 


n = 1 

s 2 c 

= B^ A sin tot + B^ A cos cot sin tot + • • • 


= — g— + oscillatory terms (H-5) 

All terms except the constant term have zero average values. Accordingly, 
if the product is time averaged over an integer number of cycles at frequency 
to, the result is: 
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(H-6) 


^ B®A 

y A sin cot = — ^ — 

s 

where the bar denotes time averaging. Solving for gives 

B® = -| y A sin cot (H-7) 

If the output y is multiplied by A cos wt (the input phase shifted 90°) and 
time averaged, a similar argument can be used to show that 

Bj = yAsinwt (H-8) 

The conceptual operation of the Fourier Analyzer is shown in block dia- 
gram form in Fig. H. 1-1. Since the required 90° phase shift is achieved 
using an integrator, the output of the quadrature component multiplier is 
the desired quantity divided by «. 


H. 2 FOURIER ANALYZER MECHANIZATION 

The analog computer mechanization diagram for the Fourier 
Analyzer is shown in Fig. H. 2-1. The analyzer is scaled in such a way 
that if A = 1 volt, the in-phase output component is one volt when the gyro 
loop is producing unity closed loop gain with no phase lag. Since the input 
scaling is 1 volt = 0. 1 rad/sec, the amplitude of the in-phase component 
B® in rad /sec can be computed: 


_ 1 /In-Phase 
' B 1 'rad/sec TO \ 


Output Component^ 


volts 


(H-9) 
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The Quadrature Correction Divisor is used to compensate for the 1/co 
factor introduced when the input signal is phase shifted. Given the Quadra- 
ture Output Component and the Quadrature Correction Divisor, the true 
quadrature component, B^, is computed. 

1 Quadrature Output Component /„ 

10A Quadrature Correction Divisor ' " ' 


^lVad/sec 



Figure H. 1- 1 Conceptual Operation of the 

Fourier Analyzer 


The mechanization is designed to prevent a constant input to the 
gyro from causing a ramp into the Quadrature multiplier. 
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Figure H„ 2-1 Analog Computer Mechanization 

of the Fourier Analyzer 


The low pass filters used to attenuate high frequency terms 
from the two multiplier outputs have a bandwidth of 0. 01 rad/sec. Con- 
sequently, the rise time of the filters in response to constant inputs is 
five minutes (three time constants). Care must be exercised to wait the 
required time period before observing the In-Phase and Quadrature Output 
Components. 
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APPENDIX I 


DERIVATION OF 

TRANSFER FUNCTIONS REPRESENTING 
ROTOR- TO- GIMBAL COMPLIANCE EFFECTS 


The use of a second order differential equation to describe the 
dynamic relation between float angle and applied torque in a single -degree - 
of -freedom gyro permits simple calculation of limit cycle and frequency 
response characteristics. The effects considered in such a dynamic model 
of the gyro float are only those of lumped rotor and gimbal inertia and linear 
fluid damping. It is now desired to investigate the effects of angular com- 
pliance between the rotor and gimbal. Angular flexure can generate sizable 
error torques about the float output axis by coupling with the spin angular 
momentum. Consequently, the potential exists for significant modification 
of the transfer functions used in determining gyro limit cycles and closed 
loop frequency response. 

Figure 1-1 illustrates the effects essential for considering rotor- 
to-gimbal compliance. The rotor is connected to the gimbal along the gyro 
input (IA) and output (OA) axes by structural members with stiffness and 
kg respectively. Structural damping is also assumed for motion about the 
axes. Angular rotations of the rotor and gimbal with respect to an inertially 
fixed reference frame are represented as A f and with the appropriate 
subscript to show axis of rotation. 
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H 



Figure 1-1 Rotor- to- Gimbal Compliance 


Summing moments on the rotor about the gyro input axis: 


I.. A = - k, 
11 r. 1 
r 1 


(a -a \+ha - cl/a -A \ (1-1) 

1 r i V r o H r i h) 


Summing moments about the gyro output axis : 




- d 2 ( A , 



(1-2) 


Assuming that the gyro case and gimbal form an essentially rigid structure 
about the gyro input axis, Ag. can be represented in Laplace Transform 
notation as the integral of the angular rate of the gyro about its input axis, 



(1-3) 
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Similarly, Ag Q can be represented in terms of the angular rate about the 
output axis and the angle between the gyro case and gimbal, a Q : 



(1-4) 


Assuming equal rotor inertias and stiffness and damping coefficients about 
the input and output axes, Eqs. (1-1) and (1-2) can be written in terms of the 
Laplace Transform operator, s, as a vector-matrix equation: 


o 

Is + ds + k 
r 

-Hs 


A 

r. 

l 


1 

+ 

441 m 

•H 

3 

l 

Hs 

2 

Is + ds + k 
r 


i 

o 

r-< 


A (k + ds) 
L J 


The rotor exerts a moment on the gimbal about the latter's output axis 
according to 


Rotor moment on gimbal = 




Using Eq. (1-5): 


_-H(k 2 + 2 kds + d 2 s 2 ) 


Rotor moment on gimbal = — k s — 

(s i + ds + If +H s 


2_2 w i 


s 2 (k + ds) (s 2 I r 2 + sdl r + kl r + H 2 ) 

(s 2 I + ds + k) 2 + H 2 s 2 g o 


(1-6) 
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For limit cycle and closed loop response analyses the variable Ag^ is 
represented by its component a 0 . A signal flow diagram showing the 
compliance effects as expressed by Eq. (1-6) is given in Fig. 1-2. The 
block marked ” Float Dynamics 11 is the transfer function between the net 
torque impressed on the gimbal and the gimbal angle, a Q . 


Float Dynamics 


1/C 

S(sTf + 1) 


where = Iqq /C. It should be noted that this differs from the definition 
of usually used. When the gyro is modeled by simple dynamics, the 
rotor and gimbal are considered to be rigidly attached and the time con- 
stant is computed using the combined moment of inertia of the two com- 
ponents. When compliance is considered, the gimbal and rotor are treated 
separately and the gimbal dynamics depend only on I 00 , the output axis 
inertia of the gimbal alone. 


Define: 



p - 2 Ve; 


(1-7) 


where v n and p are the natural frequency and damping ratio, respectively, 
of the damped second order system assumed about each axis in Fig. 1-1. 
The block diagram of Fig. 1-2 can now be reduced to a more compact form, 
shown in Fig. 1-3. The transfer function of the float dynamics, «o/®®net» 
is given by: 
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Torquer 


Figure 1-2 Signal Flow Diagram of Gyro Loop 
Including Compliance Effects 





Figure 1-3 Signal Flow Diagram Including 

Compliance Effects (Compact Form) 
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a o 

^net 


(s) 


j s 4 + 4pi' n s 3 + [2^ (2p 2 + 1 ) + w| ]s 2 + 4py 3 s + v 4 

g" a " o " _ 

g s(s + C 4 s + CgS + C 3 s +Cj + Cq) 


Iqo 


(1-8) 


where 

C 4 = 4pv n + (d/lo 0g ) + (l/T f ) 

C 3 = 2i?(2p 2 + l)+co 2 + (4pv n /r f )+(k/l 00g )+(2py n d/l 00g ) 

C 3 - 4pi^+ (tOg/ij)+ 2v n (2p + l)/lf + (2p4 1 k/l 00g )+ (d/loOg)( t/ n + 

Ci = v n + (4pu n /r f ) + (k/l 00 g)( a5 g + 

c o = (4/ r f) 

The transfer function which relates the input angular rate, Wj, to the 
moment exerted on the gimbal by the rotor, Mj_, is 

Mj -H(4p 2 i/ 2 s 2 + 4py 3 s + v 4 ) 

— -(s) = —TL X 9 9 3— 5 5 37 (I“8) 

^ s+4py n s +[ 21 ^( 1 + 2p ) + Wg]s +4py n s + i^ n 
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APPENDIX J 


RECOVERING ANGULAR RATE INFORMATION 
FROM THE OUTPUTS OF PULSE TORQUED GYROS 


The presence of a strapdown inertial navigator in a vehicle 
suggests that the same body-fixed instruments be used to provide motion 
indications for an autopilot,* Of particular interest is the vehicle angular 
rate resolved into body axes. Unfortunately, accurate strapdown gyros 
usually provide information in the form of incremental angular rotation 
about the sensor input axis and some data processing is necessary to re- 
cover an indication of angular rates. This Appendix outlines several such 
techniques. 


J. 1 ASSUMPTION OF PIECEWISE CONSTANT ANGULAR RATES 

If the body axis angular rates are varying slowly compared to 
the output frequency of the gyros, a simple scheme for recovering angular 
rate information may suffice. This involves collecting gyro output pulses 
over a given time period and dividing the net indicated angular rotation of 
the sensor by the interval during which the pulses are collected. The 

A 

estimate of angular rate about the gyro input axis, wj, is given by 


*In addition to the desire for body axis angular rate information for use in 
a vehicle autopilot, some accuracy advantage may be gained if angular 
rate information is available for direction cosine matrix calculations. 
However, this consideration is beyond the scope of the present report. 
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(J. 1-1) 


rA0 
sT 

where A0 and T are the gyro output pulse weight and clock frequency and 
r and s are integer numbers. The product rA0 is the angular increment 
measured during the period sT. In choosing the interval for calculating 

A 

coj, the product sT should be small enough to permit the approximation of 
the input angular rate as a constant. 

A 

Errors in the estimate of of Eq. (J. 1-1) arise from time- 
varying input angular rates and from quantization of the gyro output. The 
quantization error in the angular increment r A0 can be viewed as a ran- 
dom variable uniformly distributed over the range (-A0/2, A0/2) with a 
mean square value of A0 /12. The mean square error in o)j is: 

^ 1 A0 ^ 

Mean Square Error in = -rrs (=r) ; 1-isT (J. 1-2) 

Equation (J. 1-2) illustrates some factors in the accurate estimation of angu- 
lar rate from the output of pulse torqued gyros. Consistent with the assump- 
tion that Wj is constant, the interval sT should be made as large as possible. 
Also, finer gyro quantization gives smaller errors in Depending on the 
particular vehicle environment and autopilot this scheme may be perfectly 
acceptable. 

Gain and phase relations for this rate recovery technique based 
on a linear gyro model and the assumption that s >> 1, are derived in 
Ref. 8. The gain relation, expressed as a function of normalized frequency 
is illustrated in Fig. J. 1-1. The phase lag is simply -2 tt c*/co s , where 
u) s is the frequency at which new values of angular rate are computed: 
co s = 2ir/sT. From Fig. J. 1-1 it can be seen that this approach to angular 
rate recovery is satisfactory when new angular rate estimates are calculated 
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at a frequency which is an order of magnitude higher than that of the signal 
being recovered. If a limit cycle of known frequency, is present in the 
gyro output, proper adjustment of the frequency co g will remove this un- 
wanted oscillation from the angular rate estimate. For example, if 
to g = Wj i no angular rate will be indicated as a consequence of the gyro 
limit cycle after the first calculation is performed. If it is necessary to 
compute angular rate information at a higher frequency to measure trouble- 
some vehicle oscillations, some of the limit cycle will be evident and 
further filtering of the angular rate signal may be required. 



Frequency of Angular Rate Signal, Wj 
(2ir/sT) 


Figure J. 1-1 Gain Characteristics for Integrating Gyro 

and Rate Recovery Techniques Expressed 
in Eq. (J.l-1). Source: Ref. 8. 
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J.2 ASSUMPTION OF PIECEWISE 

CONSTANT ANGULAR ACCELERATION 

If the errors produced in the technique discussed in Section J„ 1 
are unsatisfactory, angular rate information can be made more sophisti- 
cated by including the possibility of angular acceleration. In this case the 
derivative of is assumed to be piecewise constant. Over any interval r 
the gyro output, rA0 is assumed to describe an angular rotation generated 
according to: 

rA0(t) = / co|(t) dt = a^t + a 2 t^ (J.2-1) 

o 

and the coefficients a^ and a2 are to be recovered from the gyro outputs. 
Two consecutive measurements of angular rotation, r n _jA0 and r n A0 , 
are required to estimate angular rate at t = nr. The coefficients aj and 
a2 are assumed constant over the time interval of length 2r. Using the 
equations: 


/'(n-l)r A A 

r R jA 0 = / Wi(t) dt = raj + 

(n-2)T 


r nT A n A 

(r 1 + r )A0 = J Wi(t)dt = 2Tai + 4r z a 2 
(n-2)r 


(J. 2-2) 


estimates of aj and a 2 at t = (n-2)t are found: 


»1 = < 3r n-l - r n) 
*2 = < r n - r n-l> 


(J. 2-3) 
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The estimate of body angular rate about the gyro input axis at time nT is 
the derivative of the right side of Eq. (J. 2-1) evaluated at t = 2 T: 

A A A 

C0|(nr) = aj + 4ra2 

Using arguments similar to those applied in Section J. 1 and assuming that 
errors in the gyro output angles indicated for consecutive periods of dura- 
tion r are not correlated, the error in the estimate of ooj caused by gyro 
output quantization has a mean square value given by 

A 2 

Mean Square Error in _ 5 1 ,A0v ,, „ 

Caused by Quantization 7Z T2' r ' ' * ” 

Comparison of Eq. (J. 2-5) with Eq„ (J. 1-2) indicates that for both approaches, 

A 

the errors in ojj have similar relative dependence on the output observation 
interval r and the gyro output quantum A0. The absolute sensitivity of the 
more complex scheme under discussion here is larger but this fact should 
be balanced against the anticipation of smaller errors from incomplete 
description of the angular rate history. 

The accuracy involved in assuming the angul ar rate model im- 
plied by Eq. (J.2-1) can be determined by expressing coj as a polynomial 
with any desired number of terms. For example: 

o>|(t) = b 0 + b 1 t + b 2 t 2 . . . (J. 2-6) 

The output of a perfect gyro over an interval r would then be 
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r(n-l)T 

^n-1 = J w i(t) dt 

(n-2 )T 

bj 9 b2 o 

— bQ^ + -!£- T + -y- T + o . a 


/ nr 

«i(t) dt 

(n-2)r 


= 2bQT + 2b jT 2 


8 . ,3 
+ S b 2 T 


(J. 2-7) 


Substituting ^ n „i for r n _j A0 and for (r n _j + r n ) A8 in Eq. (J.2-4) 
yields the estimate 

a in 

o)^(nr) = bg + 2bj? + -g- b 2 T 2 + . . . (J. 2-8) 

A 

The error in using the technique discussed in this section is the differ- 
ence between Eqs. (J.2-8) and (J.2-6): 

A 

e = aj^nr) - wj(nr) 

— bQ + 2b j r + -y b 2 T ^ + t . . 

(J. 2-9) 

- [bg + SbjT + 4b2T 2 + . „ .J 

2 9 

= --g b 2 t + higher order terms in r 

The errors caused by incomplete specification of the nature of Wj are pro- 
portional to the second and higher power of r, as would be expected from 
the fact that second and higher derivatives of input angular rate are ignored 
in this formulation. Of course, higher order polynomial representations 
of can be used or the polynomial can be fitted to more data points. For 
example, see Ref. 5. 
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J. 3 OPTIMAL FILTERING 


Optimal linear filtering can be used to recover angular rate in- 
formation from the outputs of pulse rebalanced integrating gyros. Some 
work has been performed at TASC in this area but has not been documented. 
The discussion below summarizes that effort and presents some indications 
of the accuracy which can be obtained. 

In order to apply linear filter theory to the problem of recover- 
ing angular rates , motion about the gyro input axes is described as a 
random variable. In particular, angular rate is modeled as an exponen- 
tially correlated quantity. The gyro is assumed to provide a quantized 
representation of the integral of the random angular rate. Since the mea- 
surement to be processed by the filter is quantized, recent developments 
in the design of optimal linear filters (Ref. 6) are employed to account for 
this aspect of the problem. Based on the analytic determination of steady- 
state errors in the estimate of angular rate, the normalized relations 
illustrated in Fig. J. 3-1 were obtained. Caution should be exercised in 
using these curves. The relations illustrated are based on the assumption 
that process noise acts on the system in a way that causes the conditional 
errors in estimates of the state variables to have a Gaussian distribution 
immediately prior to each measurement. This requirement depends on 
the magnitude of the process noise, the angular rate correlation period 
and the time between new inputs to the filter. 
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e 2 = Steady state mean square error in C0| 

N = Covariance of the process noise driving 
the linear first order system whose 
output represents vehicle angular rate 

q = Correlation time of vehicle angular rate 


Figure J. 3- 1 Angular Rate Estimation Errors: 

Optimal Linear Filter for Quan- 
tized Measurements 
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